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NOMENCLATURE 


A, B, C 


m 


d 

dBA 
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E 
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w 

u 

V 

X, r 
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calibration constants 

axial velocity fluctuation sensitivity coefficient 

radial velocity fluctuation sensitivity coefficient 

(E, A + E, A ) 

X Vj 2 »2 

eddy vl'jc.xslty factor 

confined jet tube diameter 

(E A + E A ) 

X ^ ^2 

nozzle exit diameter 
attenuation of channel A 
attenuation of channel B 
time-mean voltage of hot-wire = e 
voltage of hot-wire = e + e' 

radial position Indicator; j^l on the centerline 
maximum radial position Indicator 
free jet parameter 

distance from nozzle exit to contraction nozzle 

hot-wire resistance 

axial velocity component = "u + u’ 

turbulent Reynolds shear stress 

radial velocity component « v + v* • 

axial and radial coordinates 

effective cooling velocity 
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1 
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single mixing length scale 
expansion angle 
yaw angle 

turbulent viscosity 

nondlmenslonalized turbulent viscosity 
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Subscripts 
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value at X = 0 
root-mean-square value 
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Superscripts 


time-mean average 
fluctuating quantity 


CHAPTER I 


INTRODUCTION 

1.1 Investigation of Combustor Flowfielda 

A renewed Interest In gas turbine con^ustor analysis has been 
initiated by new emphasis on pollutant and noise suppression and fuel 
economy. The combustor, shown in Figure 1, see Appendix B, must ful- 
fill several requirements (1) Including: 

Bum fuel efficiently. 

J. Maintain stable combustion over a wide range of operating 
conditions . 

3. Uniform distribution of temperui.t."e in product gases. 

4. Minimum size and weight. 

5 . Minimum pressure drop. 

Flowfields within such combustors are very complex due to the corner 
and central recirculation regions produced by a rapid expansion and 
strong swirl imparted to the incoming air, llius, the flow is a swirl- 
ing, recirculating, turbulent one and is, therefore, difficult to 
analyze. Aero thermo chemistry is another formidable problem facing the 
designer, and more thorough and accurate procedures can yield the de- 
sign objectives mor«; quickly and less expensively. 

The general aim of this investigation is to provide information 
on combustor phenomena so that combustion characteristics can be deter- 
mined, These characteristics either in the fom of experimental or 
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theoretical analysis are useful to designers. Theoretical investiga- 
tors can also compare the experimental data with the data from a 
numerical solution of the appropriate governing partial differential 
equations. With this general knowledge, combustor development time 
and cost will be reduced. 

The research presently being investigated involves several methods 
of analysis for both swirling and nonswlrllng flows. Tt consists of 
flow visualization, computer simulation, and time-mean and fluctuating 
velocity measurements with the emphasis on nonreacting, steady turbu- 
lent flow in axls 3 rmmetric geometries. Two simplified flowfields are 
being investigated; a free round jet issuing horizontally into a 
quiescent atmosphere and a confined jet. In the confined jet, flow 
enters through a nozzle of exit diameter d and proceeds via a sudden 
expansion Into a tube of diameter D. A swlrler may be located up- 
stream of the inlet of the tube to swirl the floi7. Also, a contraction 
nozzle may be positioned at some downstream location in the tube. 

1.2 The Problem of Turbulence Measurement 

Complex flowfield ♦•urbulerice measurements have always been a dif- 
ficult process. Several authors have discussed turbulence phenomena in 
detail and suggested various methods of turbulence measurement (2-5). 
The hot-wire anemometer is the moat broadly used instrument to obtain 
turbulence quantities. A single hot-wire in a two-dimensional flow 
can measure the streamwlse components of the time-mean velocity and 
the root-mean-square velocity fluctuation at a particular location 
in the flowfield when used at a single orientation. By using multi- 
orientations of a single hot-wire, the three components of the mean 
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and fluctuating velocity vector can be determined (6»8). The problems 
associated with this method are that extensive data reduction is re- 
quired to calculate the desired turbulence quantities and, because the 
wire has to be positioned several times, the data obtained is not 
truely simultaneous. A multi-wire probe can be employed using a single 
orientation to determine flowfleld quantities. In a two**dimenslonal 
flow, the time-mean velocities, streaniwlse and cross stream velocity 
fluctuations and their cross correlation can be determined with a two- 
wire probe. A crossed hot-wire probe was used In this study In which 
the wires were Inclined at i 45® to the flow direction and parallel to 
the probe axis. This probe should give more accurate turbulence Inten- 
sities and cross correlation than a single wire multl-orlentatlon method 
because it Is used at only one orientation and thus the data acquired Is 
nearly simultaneous. 


1.3 Previous Experimental Studies 
1.3.1 Recent Studies at Q.S.U. 

As mentioned previously, the research presently being investigated 
Involves several methods. Still (9) and movie (10) photography of 
nf ut rally-buoyant helium-filled soap bubbles and injected smoke has 
aided in the visualization of various combustor flowflelds. A five- 
hole pitot probe has been used to measure time-mean velocities at dif- 
ferent swirl strengths (9, 11). Turbulence quantities have been 
determined with a slx-orlentatlon single hot-wire technique (12) in 
swirling and non-swlrllng flowflelds. Predictions of confined swirling 
flows corresponding to the ones ntudied experimentally have been 
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achieved using a computer code that solves the appropriate partial dif- 
ferential equations (13) . These predictions (14) are now supplemented 
by predictions deduced from downstream nozzle effects (15) and realis- 
tic Inlet conditions (16). 

1,3.2 Free Jet 

Corrsln and Uberol (17) measured the power specta of velocity and 
temperature fluctuations in heated t.nd unheated jets using a crossed 
hot-wire similar to the one used in the present study. They concluded, 
from their neasurements of the shear-correlation spectrum, that the 
fine stiucture In the shear flow In a round turbulent jet is isotropic. 

Single and double hot-wires were used by Kolpln (18) to determine 
turbulence intensity and shear stress in the mixing region of a round 
jet. He found that, above a certain Reynolds nunvbei. , the jet diameter 
plays no role In the description of the turbulent field and that the 
generation of turbulence, governed by the interaction of the large ed- 
dies and the mean flow. Is a local phenomenon. Davies et al. (19) also 
conducted experiments In the mixing region in which the turbulence in- 
tensity, turbulence spectra, and sheai stress were measured. Results 
Indicated that the local intensity of the turbulence Is equal tc 0,2 
times the shear velocity, where this 'eloclty Is defined as the product 
of the local Integral length-scale of the turbulence with the local 
shear. 

Measurements in the self-preserving region of a jet v;ere performed 
by Wygnanskl and Fiedler (20) and Rodl (21) using hot-wire probes. The 
mean velocities, normal stresses and shear stress measured by Wygnanskl 
and Fiedler were also obtained by Rodl as a check. Rod! processed the 
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signals from the hot-wire anemometer using two different methods: the 

conventional method and a new method which determines the fluctuating 
velocities from the squared electrical signal. The results agree 
between the two studies when the conventional processing is used but 
hlgh'ii^ values are Indicated in Rodl's Investigation with the squared 
data processing. Unfortunately, this new method is not applicable to 
low- intensity flows, that Is, flows with less than 25 % turbulence 
Intensities. 

Bradshaw et al. (22) conducted measurements of turbulence proper- 
ties In the noise producing region of a free round jet. It was found 
that the flow in this region Is dominated by a group of large eddies, 
containing nearly a quarter of the turbulent shear stress in the quasi- 
plane region of the shear layer. 

Corrsin (23) and Sami et al. (24) performed turbulence measure- 
ments that can be compared with the results of the present study be- 
cause of the close geometric similarity and experimental conditions. 

In his experiment, Corrsin used a four hot-wire probe to measure the 
turbulence qtiantltles and a pitot tube to measure the mean velocity in 
a heated round jet. Mean axial and radial velocities, turbulence in- 
tensities and turbulent shear were determined by Sami et al. using 
single and crossed hot-wire probes in the flow-establishment region 
of an air jet Issuing into still air. The data of Sami et al. (24) 

Is presented because of the similarity in measurement techniques. 

1.3.3 Confined Jet 

.'low separation behind an orifice in an electrically heated tube 
was Investigated by Krall and Sparrow (25) . The peak Nusselt ninnber 
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was presumed to be at the reattachment point. They concluded that the 
local heat-transfer coefficients In the separated, reattached, redevelop- 
ment regions are several times as large as those for a fully developed 
flow. Also, they found that the reattachment length was unaffected by 
Reynolds number In the turbulent range. Uils result was verified by 
Phaneuf and Netzer (26) who measured mean axial velocity and wall static 
pressure behind abrupt pipe steps. It was concluded that the reattach- 
ment zone spreads out with increasing Reynolds number and Increasing 
step height, provided the Reynolds number is In the appropriate range. 

Bade and Roschke (27, 28) measured the reattachn^nt length down- 
stream of a pipe step by dye visualization In water flows, A conical 
contraction section was placed just upstream of a sudden expansion to 
reduce the boundary layer thickness. Dye was Injected through holes 
positioned along the larger diameter tube at Intervals of one step 
height. Reattachment locations were then observed by Inspecting the 
direction of dye flow. Visual observation of small buoyant particles 
in a fluid in the redevelopment region was accomplished by Feuersteln 
et al. (29). The length of reattachment found by Feuersteln Is slight- 
ly different than predicted by Back and Roschke probably due to dis- 
similar Inlet profiles. 

Pitot probe and hot-wire measurements in various sudden expansion 
flowflelds were made by Ha Mlnh and Chassalng (30). The time-mean 
axial veloclt'' was determined with both Instruments and good agree- 
ment was exhibited. A rotating. Inclining, single hot-wire technique 
was employed ‘;o measure the Reynolds stresses. One of their expansion 
rctlos matched the one used In the present study but the axial loca- 
tions of their measurements were not compatable. 
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Hoon and Rudinger (31) and Johnson and Bennett (32) developed 
laser-Doppler Systems to measure mean velocities and recirculation re- 
gion geometries, but the expansion ratios differed from the one used In 
the present experiments. 

Because of the differences In geometry, comparisons between the 
results of the cited experiments and the present experimental study are 
not possible. However, the experiments of Chaturvedl (33) are com- 
parable. Mean velocity and pressure downstream of a sudden expansion 
of a diameter ratio of 2.0 were measured with a stagnation tube and 
pitot probe. A single and crossed hot-wire were used to measure the 
turbulence intensities and turbulent shear stress. Mean velocity was 
also measured with the single hot-wire. Also comparable are six- 
orlentatlon single hot-wire measurements recently obtained by Jackson 
and Lllley (34) In the same test facility as in the present study. 

1.4 The Present Contribution 

To further develop combustor flowfleld prediction techniques, in- 
cluding turbulence modeling, there is a strong need to obtain experi- 
mental estimates of the mean and turbulent flow quantities and the 
turbulent viscosity. In the present study, a crossed hot-wire single 
orientation technique Is employed to measure mean velocity, turbulence 
intensities and Reynolds stress in a round free jet and in a combustor 
simulation confined jet flowfleld. The free jet measurements provide 
an independent check on the data acquisition and interpretation techni- 
ques. Measurements in the confined jet are carried out for nonswirling 
flow with a sudden expansion. 

Chapter II describes the experimental facilities and hot-wire and 
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raw data acquisition instrumentation. This chapter also explains the 
concept of fluid flow at abrupt expansions. Experimental techniques, 
data reduction scheme and turbulent visc:?9lty formulation are f»lven in 
Chapter III. The computer program utilized in the data reduction and 
user's guide are contained in Appendixes D and C respectively. Mea- 
surement procedures and results for free jet and confined sudden ex- 
pansion flows are discussed in Chapter IV, Tables of output quantities 
are contained in Appendix A and Appendix B consists of the figures. 

Some of the turbulence quantities are compared with measurements done 
by Sami et al. (24) in the free jet and Chaturvedl (33) and Jackson 
and Lllley (34) in the confined jet facilities. Chapter V concludes 
by summarizing the results of the present study and suggests further 
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CHAPTER II 


EXPERIMENTAL FACILITIES 
2,1 Idealized Flot<rflelds 

A round free jet issuing into a still ambient fluid, shown in 
Figure 2, has been the subject of many experimental investigations. The 
flow may be separated into three regions, each having a distinct set 
of characteristics. Region I contains a potential core of uniform mean 
velocity bounded by a shear layer. This region starts at the nozzle 
exit and extends to four or five nozzle diameters downstream. In 
region II, the velocity on the jet centerline decreases with increasing 
distance. Region II Is also characterized by the lack of similarity in 
the Tuean velocity distribution. After about eight nozzle diameters, 
region III stares and the jet flow becomes fully developed in which the 
mean velocity flow patterns in consecutive sections become similar. 

Confined flow at abrupt expansions, as in an actual gas turbine 
combustor. Is a typical example of the problem of fluid flow separation 
from the boundary. A region of reverse flow is caused by this separa- 
tion at the sudden expansion. This region, called the comer recir- 
culation zone (CRZ), is associated with a large pressure drop which 
adversely effects the performance of a can-type combustion chamber. 

The turbulence levels are high on the shear layer between the two 
regions which Indicates that, at this position, good combustion would 
occur. Flows In actual reacting can combustors have an area of high 
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heat transfer at the location of the reattachment of the dividing 
streamline which could lead to rapid deterioration of the walls of the 
combustor. A schematic of the test section geometry and expected re- 
circulation zones are shown in Figure 3. The central torroldal recir- 
culation zone (CTRZ) is present only with swirling flow. 

2.2 Free Jet 

The free Jet facility, shown in Figure A, consists of a seamless 
contoured nozzle fed by a thermally stabilized compressed air generator 
which delivers the desired flow rate through a small pressure regulator 
and a Fisher and Porter Model 10A1735A rotameter. The nozzle has a 
throat diameter of 3.4 cm and an effective flow management section just 
upstream. 

It is necessary to calibrate the crossed hot-wires in a flow of 
knomi characteristics. The potential core of the free Jet was utilized 
for this calibration. The sensitivity to angulation (see Section 3.1) 
was obtained by the use of a rotary table on the jet centerline to 
rotate the probe inside the potential core from -10“ to +10“ to the 
free stream flow direction. An "L" shaped probe bracket holds the 
probe while it is being calibrated. 

2.3 Confined Jet 

The Oklahoma State University Confined Jet Test Facility is used in 
this experimental study. It 1s a simulation of a typical axisymmetrlc 
combustion chaidier of a gas turbine engine. Figure 5 shows a schematic 
of the overall facility. Ambient air enters the low-speed wind tunnel 
through a foam rubber air filter. The air then flows through an axial 
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flow fan driven by a 5 h.p, varldrlve motor which allows the flow rate 
to be varied for different test conditions. Next, the flow passes into 
a flow conditioning section which reduces the intensity of the turbu- 
lence caused by the driving fan. This turbulence management section 
consists of a perforated aluminum plate, followed by a fine mesh screen, 
a 12.7 cm length of packed straws, and five more fine mesh screens. 

The flow then enters a contoured nozzle designed by the method of Morel 
(35) to produce a minimum adverse pressure gradient on the boundary 
layer to avoid local separation and flow unsteadiness. The ratio of 
the turbulence management cross section to that of the nozzle throat, 
by area. Is approximately 22.5. The nozzle throat diameter d Is 15 cm. 

Finally, the flow enters the test section, which Is an idealized 
combustion changer model. It Is composed of a swlrler (optional), an 
expansion block and a long plexiglass tube. The expansion block Is a 
30 cm diameter disk of wood that has a 15 cm diameter hole centered 
on its axis and is attached after the swlrler. In this study, the 
swlrler Is not used. Thus, the expansion block is mounted directly 
onto the throat of the wind tunnel nozzle. The flow expands Into a 
plexiglass tube of diameter D of 30 cm, thus giving a diameter ex- 
pansion ratio (D/d) of 2, The test section Is constructed of plexi- 
glass for flow visualization and ease of location of measuring probes. 
Holes are located at multiples of x/D <= 0.5 along the length of the 
tube, so that the hot-wire probe can radially traverse the airflow. 

A contraction nozzle with an area ratio of 4, shown In Figure 6, may 
be positioned downstream of the expansion block to simulate the con- 
fining wall of a real gas turbine coid>U8tor. The test section does 
not have film cooling holes or dilution air holes, typical in a real 
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combustor. The ample size of the test section used in this study pro- 
vides good probe resolution for the hot-wire measurement a . 

2.4 Hot-wire Instrumentation and A/D Conversion 

As mentioned earlier, a crossed hot-wire probe, DISA Type 5SP63, 
is used in this experimental study in which the wires are inclined at 
1 45° to the flow direction and parallel to the probe axis. The probe 
orientation in the flowfield is illustrated in Figure 7. Each hot- 
wire has a separate constant temperature anemometer for complex flow 
measuLement. A circuit diagram for the DISA 5SM01 main frame with a 
55M10 standard bridge anemometer is given in Figure 8. To obtain the 
turbulence intensities in the axial and radial directions, it is neces- 
sary to Instantaneously add and subtract the anemometer voltage outputs. 
The adder and subtractor circuits were constructed of commercially 
available 741 op amps and are shown in Figure 9. The multiplication of 
the voltages by a Salcor Model SAI 43 correlation and probability ana- 
lyzer will give the shear stress. Because the fluctuating signals are 
small, each signal is amplified by a Hewlett Packard Amplifier with a 
20 dB gain. The frequency response of the hot-wires and the associated 
electronics is approximately 40 kHz based on square wave response teats. 
ITie use of llnearlzers was judged unnecessary from experience in the 
use of crossed hot-wire anemometry gained previously at Oklahoma State 
University by researchers such as Iltorrison (36) and Swearingen (37). 

To Increase the accuracy of the measurements, a microprocessor 
with 12-blt analog to digital (A/D) conversion is employed to acquire 
the data. A schematic of the hot-wire and data acquisition equipment 
layout la shown in Figure 10. Because the A/D converter, an Octogon 
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Sys-2, has a maximum Input voltage of 2.5 V, the Incoming mean voltages 
are stepped down by a factor of four. The nos signals from the adder 
and subtracter are changed from AC to DC signals so that the A/D conver- 
ter can read them. Accuracy with 12-blt conversion is: 

' ^ ' X 2,5 V X 4 = 0.0024 V for mean voltages 

2*2_i 


"" ■ X 2.5 V - 0.006 V for fluctuation voltages 

2^2_i 

The 8anq>llng rate Is 100 per second if the microprocessor with A/D con- 
version is just reading the Input signals. In this experimental study, 
however, the microprocessor calculates an average of 625 voltages for 
each Input signal and thus the sampling is impared to the point of 
reading 12 voltages per second. 

A traversing mechanism shown in Figure 11, supports the hot-wire 
in the confined jet test facility. The hot-wire probe is Inserted 
Into the flow through a rotary vernier and a base that mounts on the 
outside of the plexiglass tube. With the rotary vernier, it is possi- 
ble for the probe to be traversed to any radial position at selected 
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CHAPTER III 


EXPERIMENTAL TECHNIQUES AND 
DATA REDUCTION 


3.1 Experimental Techniques 


The hot wires are assumed to obey an extension of King's Law. lliat 


is : 


1 

E = A + BZ"^ + CZ 


where A, B and C are constants determined a priori by a calibration 
experiment and Z is the effective cooling velocity. Thus for each wire: 



*1 * + <=i"i 


^2 ^ S ^2 


( 1 ) 


The calibration is accomplished prior to every experiment by plac- 
ing the crossed hot-wire on the free jet centerline within the potential 
core region. The flow rate through the nozzle and thus the exit velo- 
city is varied by an upstream control valve and determined by a rota- 
meter. For each rotameter setting, there is a corresponding axial 
velocity. Thus the voltage from each of the wires is recorded for 
each rotameter setting and calibration curves of voltage, E, versus 
velocity, "u, can then be drawn. These are presented in Figure 12. 

If Equation (1) is Inverted, the effective cooling velocity on 
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each wire is obtained, 

= [(-B^ + - 4C^(A^ - 

Z 2 “ C(-B2 


Therefore, with a p-'Ban voltage from wire 1 or 2, the corresponding 
velocity can be found. Using the geometry shown in Figure 7, the mean 
axial and radial velocities are given respectively by: 


V 



Tlieae velocities are nondlmenslonallzed with respect to the nozzle exit 
velocity, u^. 

The instantaneous velccity fluctuations measured from an inclined 
hot-wire can be defined by the '"allowing expression derived by Corrsin 
and Uberoi (17) ; 


E 


A u* ± A V* 
m ~ V ^ 

u u 


( 2 ) 


where A^ is the sensitivity coefficient for axial velocity fluctuation 
and A^ is the radial velocity fluctuation sensitivity coefficient. 
These coefficients are evaluated in the following manner: 


A 

m 


a InE 
a Inu 


R , (b constant 
w 


A 


V 


a InE 

a^) 


R , u constant 
w 
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If a crossed hot-wire Is used at the orientation shown in Figure 7, then 
equation (2) written for each wire becomes: 

e ’ - E A ~ + E, A ^ 

1 1 "lu 1 ''lu 


e ' - E A a~ - ^ k ^ 

2 ^ "2u 2 Vj- 


Output from each wire can be Instantaneously added and substracted to 
obtain: 


u» . . v' 


(e/ + e.') = (E,A + E A ) — + (E,A - E A ) — 

1 2 1 2 m2 - 1 Vj_ 2 V 2 - 


(e • - e •) = (E A - E A ) (E A + E A ) ~ 

^ ^ ^ "1 2 m2 - 1 2 V2 - 


The crossed hot-wires are matched so that the sensitivity coefficients 

A and A are approximately the same for both wires and the overheat 
m V 

ratios of the wires are set so that E^^ = E^. These simplifications 
give: 



*'a'> ■ 

U 


where D = 


m 


(E.A 
1 m- 


ni/i 


and 


(ei- - - B ^ 

U 


where B, = (E,A + E A ) 

•Jt 1 2 V2 

Thus, solving for the turbulent velocities 
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u 




m 


U (e^- - 


which may be nondlmenslonallzed via 


u'/u “ (u'/tT) (u/u ) 
o o 

v'/u = (v'/u) (“ir/t ) 
o o 


To obtain the 


A and A need to 

m V 


appropriate D 

m 

be evaluated. 


or B,, the sensitivity coefficients 

<P 

That is; 


9 InE u 9 e 

t = — — — a — 

™ 9 Inu E 9u 


9 InE ^ jL 3E 
3(() E 


To determine A and A , 9E /9u and 
”l “2 ^ 

derivative of E with respect to u of 


BE^/Su are calculated from the 
Equation (1). Thus, 


9 E jv *5 ^ 

i ~ [(A. + B u + C u) 

9u 9u ^ ^ 

2C^n + 

J J- ig 

4(Aj^u + B^u + Cj^u ) 
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Slmllarlly for wire 2: 



h 

2C^u + B 2 

- — 

+ B2~ + C^u )'■ 


Therefore, an A and A can be determined for a certain number of 

®2 

velocities and with these, a graph of D against velocity can be drawn. 

ra 

The relationship between D and velocity is a linear one and is shown 

m 

in Figure 13. Thus any D la knovm within the range of velocities 

m 

encountered. 


The sensitivity to angulation coefficients A and A are obtain- 

^1 "^2 

ed by positioning the probe in the X-Z plane of the round free Jet and 
rotating the probe in 1° increments between +10° and -10° yaw angle to 
the free stream flow direction, shown in Figure 14, Both mean voltages 
are recorded for each yaw angle <!’ and a linear curve fit Is then ap- 
plied to each voltage versus Figure 15 contains examples of these 
measurements. Because the wires are matched to a high degree, the 
curves Intersect at or very near the zero degree point. The slopes of 
these plots are 3E^/d<j> and dE^/d<f>» This procedure is repeated at the 

same velocities at which the D ’s are calculated. Thus, the sensitivity 

m 

coefficients A and A can be determined for each angulation test 

Vi V2 

velocity. A graph of against velocity can then be plotted and again 
can be approximated by a straight line as illustrated in Figure 13. 

A cross correlator can be used to determine the major shear com- 
ponent of the Reynolds stress tensor in a turbulent shear flow by 
multiplication of the voltages proportional to u*/u and vVtl. Thus: 
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Correlator Output 

Bx • 

(p m 



From the Honeywell SAI-43A correlator manual; 

CrreLtor Output - 2^2^ 

1Z5 

where X Is the voltage correlation output, 

dBA Is the attenuation of channel A, and 

dBB is the attenuation of channel B. 

2 

A factor of (40) is Introduced because both hot-wire fluctuating signals 
are amplified by a gain of 40. Therefore: 

uW^ , X(VoltB)riO^^°^ u 

2 125 . 1600 • Bx • D 2 

u 9 m u 

o 0 

A computer code was designed to accomplish the reduction of the 
voltage data from the hot-wires In the technique just discussed. Appen- 
dix C contains the User's Guide to the program with sample output listed 
In Appendix D. 


3,2 Turbulent Viscosity Formulation 

While the molecular viscosity Is a real property of the fluid, exist- 
ing whether the fluid Is in motion or at rest, the eddy viscosity requires 
some flow of fluid to become effective and la thus not a property of the 
fluid but a parameter of the fluid motion. This parameter describes the 
behavior of the turbulent stresses in terms of the mean-velocity gradi- 
ents Implying the assumption that the turbulence transport Is of the 
gradient- typ e . 
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The first-order closure (seml-emplrical) methods have been applied 
extensively to the calculation of the eddy viscosity. The first techni- 
que, according to the Bonsslnesq (38) concept, asserts that the eddy 
viscosity Is a scalar quantity (much larger than the laminar viscosity) 
that relates turbulent shear stress to the time-mean velocity gradients 
via; 


« -u^v** / (3u/3r + 3v/3x) (3) 

and generally varies throughout the flowfleld. 

Prandtl (39) proposed that the eddy viscosity could be specified in 
terms of local parameters; a length scale and the time-mean velocity 
gradient. He suggested; 


V 

t 


C 

y 



where is a constant and 

1 is the mixing length. 

Both methods require that the flow is not far from equilibrium and 
that the turbulence Reynolds number is high. The near equilibrium as- 
sumption has been shown to be approximately correct for a mixing layer 
(40) and the Reynolds number at which experiments are conducted in this 
present study are sufficiently large. 

In the present study, measured turbulent shear stress and time- 
mean velocity gradients are used to deduce values of eddy viscosity via 
the Bousslnesq concept. Because the magnitude of the mean radial 
velocity V is small compared to the mean axial velocity u^ Equation (3) 


becomes: 



^ 
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° -u'v* / 3Ti/3r 

where u'v' Is the turbulent Reynolds stress and ^u/^r is evaluated by 
finite difference techniques. Thus. 

3u/3r « - u^) 

3u/3r » — (Uj^i - Uj_^) 

3u/3r ” -S7 " “j-1^ 

From Inverse analysis of experimental time->mean measurements. Hlnze 
(2) states that the tuibulent viscosity in a rpund free jet is approxi- 
mately constant In the central part of the fully-developed region and 
given by the relation: 


for j = 1 


for j >j>l 
max 


for j = j 


max 


V “ 0.00196 K u d 
t o 

where K is a parameter between 3.4 and 6.39. The latter value is used 
in agreement with earlier measurements (41). In the free jet, a non- 
dlmenslonallzed turbulent viscosity v^* given by: 

\> * o \) / (6.39 u d> 
t t o 

should approach the value of 0.00196 in the fully-developed region. 

Later results :^11 confirm this. In the confined jet. the simple normali- 
zation 

V * cj V /(u d) 
t to' 


is used to be consistent with on-going theoretical studies (42). Its 



CHAPTER IV 


RESULTS AND DISCUSSION 
4.1 Measurement Procedures 

Nonewlrllng nonreacting free and confined jet flows, with coordinate 
systems shown In Figure 16, are Investigated by the use of a crossed 
hot-wire to measure mean velocities u and v, fluctuating velocities u’ 
and v', and Reynolds stress u'v' . The proceeding quantities are normali- 
zed with the nozzle uniform axial exit velocity u^, determined with the 
probe from an independent measurement at the nozzle exit face. Radial 
distribution of turbulent scoslty Is also deduced. The Reynolds num- 
bers of the flows under Investigation are high enough to ensure that the 
results are In the Reynolds number independent regime. 

Prior to any measurements, the crossed hot-wire probes were cali- 
brated in the free jet facility as described in Section 3.1. The 
voltage divider box and A/D converter were also calibrated by use of 
a voltage standard. This was to check that no 'drifting* of the elec- 
tronics had occurred. 

Measurements were made in a free jet at axial locations of x/d of 
1, 6, 10, and 20. Because published data exits for the free jet using 
hot-wire anemometry. It was used to validate the experimental procedure 
and data reduction discussed in Chapter III since little published data 
is available for confined jet flowflelds using the present method and 
test section geometry. The results are compared with those of Sami et 
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m ■ 
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al, (24). These results in tabular form are contained in Table I. 
Appendix A. 

Measurements were also made in a model of an axlsymmetrlc can com- 
bustor (confined jet test section) with an expansion ratio D/d=^2. an 
expansion angle ct of 90° and axial location Increments of 0.5 chamber 
diameters. The axial location of zero is actually 1 mm downstream of 
the enlargement face so as to prevent the hot-wire probe from damage. 
Comparable are measurements by Chaturvedl (33) and Jackson and Lilley 
(34). A single hot-wire multi-orientation technique was used by Jackson 
and Lilley (34) in the same test facility. Table II consists of flow- 
field data for the confined jet without contraction nozzle. Also, a 
confined jet flow field was Investigated with a contraction nozzle of an 
area ratio of 4 located at L/D=2. The data for this flowfield is con- 
tained in Table III. This nozzle, mentioned earlier, has a 45° slope 
facing upstream. Time-mean flowfield data for this configuration have 
been presented prevlotisly by Yoon and Lilley (43), using a five-hole 
pitot probe technique, 

4,2 Free Jet Results 

4.2,1 Mean Velocities 

Presented In Figure 17 are tjeasurements of time-mean axial velocity 
u and radial velocity v. These resiilts Indicate that the present jet 
is similar but slightly narrower than the one used by Sami et al, (24). 
The comparability of the jets also indicates that the experimental 
techniques and data reduction are sound for time-mean axial velocity 
but the time-mean radial velocity distribution shows large scatter in- 
dicating a slight error in the data reduction or experimental techniques. 
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4.2.2 Turbulence Intensities 

Figure 18 the radial distribution of axial and radial turbu- 

lence Intensities, Fair agreement with Sami et al. (24) can be seen In 

u* /u , The values measured here tend to be lower than those pre- 

nns o 

sented by Sami et al. (24) although the trends are alike. Tlils can be 

explained by the difficulty In measuring turbulent, low speed flows. 

Sami et al. did not present v' data. However, they asserted that 

rms 

u' = yr V* . This relationship is found to be approximately true 
rms rms 

in the present study. Moving further downstream, the distribution of 
the turbulence Intensities become less pointed and more flat. This 
Indicates that the turbulence levels are spreading laterally. 

4.2.3 Shear Stress 

The plots depicting the turbulent, shear stress are shown in Figure 
19. Again, fair agreement with Sami et al. (24) is to be noted. The 
shear stress is the most difficult to experimentally measure. This 
probably accounts for the slightly higher peak as compared with Sami et 
al, (24) at x/d^b. As before in the plots of the normal stresses. 
Figure 19 shows the dissipation of turbulence throughout the flowfield. 
The radial distributions of normal and shear stress also indicate that 
the present Jet is thinner than the one used by Sami et al, and that 
the nozzle Is of slightly different design causing variations in the 
stresses as compared to the present study. 

4.2.4 Turbulent Viscosity 


The radial distribution of turbulent viscosity calculated from the 
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shear eccesB and axial velocity gradient Is shovm In Figure 20. As 
can be seen, the turbulent viscosity grows from a splkey, small magni- 
tude profile to one with a larger magnitude and flat appearance. This 
flat type of distribution with a magnitude of 0.00196 Is deduced by 
Hlnze (2) when In the similarity region. The large scatter of the 
turbulent viscosity at x/d«10 and 20 Is due to the high sensitivity of 
to a change In the time-mean axial velocity, 

4i 3 Confined Jet Results 

4.3.1 Mean Velocity 

Figure 21 shows measured values of time-mean axial and radial velo- 
cities. A nearly flat axial velocity profile is seen at the entrance 
of the test section. Comparisons of the time-mean axial velocity with 
Chaturvedl (33) and Jackson and Lilley (34) show good agreement except 
in the expected region of recirculation. This Is due to two reasons: 
although a hot-wire anemometer cannot sense flow direction, Chaturvedl 
(33) and Jackson and Lilley (34) plotted some negative values of time- 
mean axial velocity by deducing direction from earlier pitot probe data 
and because of the ”L" shaped probe design (see Figure 14) the hot-wires 
are hidden from the recirculating flow by the probe supports. Because 
the time-mean radial velocity is calculated from the subtraction of two 
large qt^antlties producing a small number, the profiles shown In Figure 
21 are almost zero at all axial locations, Bruun (44) states that velo- 
cities determined by a crossed hot-wire probe have a very strong yaw 
dependence, particularly This also may account for the error in the 

time-mean radial velocity. Thus, the experimental technique is not 
adequate for these measurements. 


The time-mean velocities in the confined jet with the contraction 
nozzle at L/I^2 are presented In Figure 25. The profiles show that tlie 
results change very little from the Identical flowfleld without a con- 
traction nozzle, as also found elsewhere (34,43). 

4.3.2 Turb ulence Inteneit ies 

Inspection of the radial distributions of u’ /u and v* /u dls 

nns o rms o 

played in Figure 22 shows similar trends to those found earlier (33, 34) 
The axial turbulence intensity values are lower than previous measure- 
ments, but the radial turbulence intensity values are in good agreement. 
The axial turbulence intensity Is slightly larger than the corresponding 
radial quantity with maximum values on the shear layer between recircu- 
lating and nonrecirculating flow. Figure 22 also shows that the radial 
turbulence intensity is damped to a higher degree than the axial turbu- 
lence intensity by the presence of the test facility tube wall. 

Figure 26 presents turbulence intensities with the contraction 
nozzle and shows little change from the non-contracted counterpart. 

4.3.3 Shear Stress 

Shown in Figure 23 is the radial distribution of turbulent shear 
stress. As can be seen, Chaturvedi's (33) data indicates lower values 
of shear stress than the present study. Also, results are similar to 
Jackson and Lllley's (34) with the present measurements having smoother 
profiles. The shear layer, which is very thin at the entrance to the 
test section, grows quickly to encompass a large amount of the flow- 
field. All shear stress values are plotted as positive, although In 
the vicinity of the recirculation zone, the gradient of the time-mean 
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axial velocity in the radial direction might be positive, Infetring 
negative values of u*v*. 

Turbulent shear stress profiles of the confined Jet flowfield with 
contraction nozzle in place are presented in Figure 27. Lilce the normal 
stresses, there is insufficient change with respect to the flowfield 
without the contraction nozzle to warrant any new discussion. 

4.3.4 Turbulent Viscosity 

Figure 24 shows the radial distribution of turbulent viscosity. 

As shown, the turbulent viscosity grows in magnitude while progressing 
downstream with large peaks near the centerline (r/D =0) from 1,0 until 
2.5 diameters. Here the peak reduces because the mean velocity gradi- 
ent near the centerline has increased. Although the change in mean 
axial velocity is little at small radial and axial positions, the peak 
is not seen since the shear stress is correspondingly small, 'fhe peaks 
in the vicinity of r/D ® 0,25 at axial locations of 0.5 and 1.0 give 
an indication of the position cf the recirculation zone. The expected 
free jet turbulent viscosity of 0.00196 corresponds to 0,5 on the hori- 
zontal axis. As can be seen, the values of viscosity are lower than 
this expected value in the low shear regions and larger In high shear 
regions. It must be noted that the turbulent viscosity is the most 
difficult of all the turbulence quantities presented to ascertain. 

The radial distribution of turbulent viscosity in the confined 
jet with contraction nozzle is presented in Figure 28. As before, 
the results indicate that the contraction nozzle effects little 
change. 



CHAPTER V 


CLOSIFRE 

5.1 Suiranary and Conclusions 

A crossed hot-wire probe has been enjployed to measure the time-mean 
and fluctuating velocities and shear stress in nonswlrllng nonreacting 
free and confined jet flot^s. The turbulent viscosity is also calculated 
from these results. 

The free jet Investigation was used to validate the experimental 
technique and data reduction. The results of the free jet measurements 
show good conq>arlson with available data. 

Measurements v;ere also made in a model of an axisymmetric can com- 
bustor (confined jet test sc-’ction) with an expansion ratio D/d«2, an 
expansion angle, a of 90° and axial location increments of 0.5 diameters. 
A contraction nozzle was located at L/D=2 for some of the confined jet 
measurements. Good agreement with previous workers is to be noted. The 
confined jet investigation shows that the crossed hot-wire probe cannot 
accurately measure recirculating flow without prior direction knowledge, 
and that the experimental technique is inadequate for deduction of the 
time-mean radial velocity. The crossed hot-wire technique should give 
more accurate results of the turbulent shear stress and thus turbulent 
viscosity than a multi-orientation single-wire technique. Smoother 
profiles for turbulent shear stress are seen in the present study than 
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In the previous single-wire study. Faster development of velocitv pro- 
fileSf higher values of turbulence quantities, and larger time-mean 
radial velocities produced by the recirculation zone are the effects 
of confining the free jet. Also, time— mean and turbulence characteris- 
tics with the contraction nozzle at L/D=2 show little change from that 
of the corresponding flowfield without a contraction nozzle. This is 
probably due to the hl^ time-mean axial velocity pushing the effects 
of the contraction nozzle downstream. 

5,2 Recommendations for Further Work 

Work continuing with the crossed hot-wire should concentrate on 
deriving a new method to determine the time-mean radial velocity and 
the yaw dependence factor. Also, a new crossed hot-wire probe with a 
different geometry or rotating the probe 180 degrees in recirculation 
regions might alllevlate the flow reversal insensitivity problem. If 
the probe could be set on the main streamline at a particular location, 
all turbulence quantities with respect to the probe could be found 
assuming that the probe could be rotated. Then appropriate coordinate 
transformations could relate the quantities in the probe coordinate 
fraiiffi to quantities in the facility coordinate frame. 
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TABLE I 

FREE JET FLOWFIELD DATA 


x/d 


r/r^ 

1 

6 

10 

20 

0 

1.01100 

0.95400 

0.64800 

0.32700 

0.1 

1.01300 

0.94400 


0.32300 

0.2 

1.01400 

0.92700 

0.63600 

0.31600 

0.3 

1.01500 

0.90100 



O.A 

1.01100 

0.87600 

0.60000 

0.30600 

0.5 

1.00900 

0.83600 



0.6 

0.99800 

0.79500 

0.55700 

0.30600 

0.7 

0.96900 

0.73900 



0.8 

0.87600 

0.69400 

0.51GU0 

0.30000 

0.9 

0.72000 

0.64100 

0.49600 


1.0 

0.52300 

0.59200 

0.47100 

0.26600 

1.1 

0.36200 

. 0.52400 

0.44400 


1.2 

0.21200 

0.47000 

0.42500 

0.27300 

1.3 

0.13500 

0.43400 



1.4 

0.08100 

0.40400 

0.37000 

0.24100 

1.5 

0.05500 

0.31600 



1.6 


0.29200 

0.31900 

0.22700 

1.7 


0.30500 



1.8 


0.25200 

0.27900 

0.21300 

1.9 


0.21200 



2.0 


0.16700 

0.23600 

0.20700 

2.2 



0.20400 

0.20900 

2.4 



0. 17300 

0.17100 

2.6 



0.14600 

0.16500 

2.8 




0.14600 

3.0 




0. 13700 

(a) 

u/u 

0 





TABLE I (Ck)ntlnued) 


x/d 



1 

6 

10 

20 

0 

- 0.01000 

- 0.00700 

- 0.00100 

- 0.00500 

0.1 

- 0.00800 

- 0.00400 


0.00900 

0.2 

- 0.00300 

0.00000 

0.00500 

- 0.00400 

0.3 

- 0.00100 

- 0.00100 



0.4 

- 0.00600 

0.00800 

0.00800 

0.00600 

0.5 

0.00200 

0.01400 



0.6 

0.00700 

0.00900 

0.01700 

0.00600 

0.7 

0.01500 

0.01000 



0.8 

0.02000 

0.01200 

0.00500 

0.01000 

0.9 

0.02400 

0.02100 

0.00900 


1.0 

0.01000 

0.01700 

0.01300 

- 0.00300 

1.1 

0.00000 

0.00600 

0.00700 


1.2 

- 0.01800 

0.00900 

0.00400 

0.00100 

1.3 

- 0.01700 

0.01400 



1.4 

- 0.01400 

0.02600 

0.00800 

- 0.00600 

1.5 

- 0.01100 

0.01300 



1.6 


0.00100 

- 0.00100 

- 0.01000 

1.7 


0.00700 



1.8 


0,01400 

0.00400 

0.00000 

1.9 


0.00500 



2.0 


- 0.00500 

0.00400 

- 0.00300 

2.2 



0.00200 

0.00000 

2.4 



0.00300 

0.00000 

2.6 



- 0.00500 

- 0.00300 

2.8 




- 0.00100 

3.0 




0.00000 


TABLE I (Continued) 


x/d 



1 

6 

10 

20 

0 

0.02000 

0.09100 

0.11600 

0.06400 

0.1 

0.02100 

0.09000 


0.06400 

0.2 

0.02400 

0.09500 

0.11700 

0.06300 

0.3 

0.02700 

0.10000 



0.4 

0.03300 

0.11600 

0.11600 

0.06200 

0.5 

0.04200 

0.12000 



0.6 

0.05700 

0.12600 

0.11400 

0.06200 

0.7 

0.08100 

0.12900 



0.8 

0.11400 

0. 13300 

0.10800 

0.06300 

0.9 

0.13800 

0.13000 

0.10800 


1.0 

0.12600 

0.12600 

0.10600 

0.05700 

1.1 

0.09800 

0.12100 

0.10100 


1.2 

0.06000 

0.11400 

0.10000 

0.05700 

1.3 

0.03800 

0.10700 



1.4 

0.01900 

0.10100 

0.09000 

0.05200 

1.5 

0.00900 

0.08500 



1.6 


0.07900 

0.08100 

0.05100 

1.7 


0.08000 



1.8 


0.06900 

0.07200 

0.04900 

1.9 


0.06100 



2.0 


0.04700 

0.06200 

0.04500 

2.2 



0.05400 

0.04500 

2.4 



0.04600 

0.03900 

2.6 



0.03900 

0.03800 

2.8 




0.03300 

3.0 




0.03000 

( c ) 

u * /u 

mis o 
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TABLE X 

( Ckmtlnued ) 



1 

6 

x/d 

10 

20 

0 

0.00800 

0.06500 

0.08800 

0.05100 

0.1 

0,01000 

0.06500 


0.04900 

0.2 

0.01100 

0.06800 

0.09000 

0.05000 

0.3 

0.01600 

0.07200 



0.4 

0.02100 

0.07700 

0.08700 

O . O^rtJuO 

0.5 

0.03000 

0.08200 



0.6 

0.04400 

0.08700 

0.08700 

0.04700 

0.7 

0.07200 

0.08900 



0.8 

0.10500 

0.09200 

0.08200 

0.04800 

0.9 

0.12800 

0.09200 

0.08300 


1.0 

0.11900 

0,09100 

0.08000 

0.04500 

1.1 

0.09100 

0.08900 

0.07900 


1.2 

0.05200 

0.08400 

0.07600 

0.04500 

1.3 

0.02900 

0.08100 



1.4 

0.01300 

0.07600 

0.06900 

0.04100 

1.5 

0.00700 

0.06200 



1.6 


0.05900 

0.06100 

0.03900 

1.7 


0.06000 



1.8 


0.05100 

0.05500 

0.03600 

1.9 


0.04200 



2.0 


0.03300 

0.04700 

0.03500 

2.2 



0.04000 

0.03500 

2.4 



0.03300 

0.02900 

2.6 



0.02700 

0.02800 

2.8 




0.02500 

3.0 




C . 02300 

( d ) 

V* /u 

rms 0 
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TABLE I (Continued) 


x/d 



1 

6 

10 

20 

0 

0.00000 

0.00000 

0.00097 

0,00115 

0.1 

0.00000 

0.00000 


0.00092 

0.2 

0.00000 

0.00000 

0.00418 

0.00158 

0.3 

0.00000 

0.00338 



0.4 

0.00000 

0.00539 

0.00590 

0.00166 

0.5 

0.00000 

0.00852 



0.6 

0.00000 

0.01071 

0.00689 

0.00143 

0.7 

0.00093 

0.01143 



0.8 

0.01115 

0.01163 

0.00696 

0.00207 

0,9 

0.01296 

0.01067 

0.00692 


1.0 

0.00896 

0.00962 

0,00671 

0.00229 

1.1 

0.00477 

0.00831 

0.00615 


1.2 

0.00155 

0.00700 

0.00598 

0.00190 

1.3 

0.00042 

0.00594 



1.4 

0.00003 

0.00518 

0.00487 

0.00194 

1.5 

0.00000 

0.00376 



1.6 


0.00326 

0.00383 

0.00191 

1.7 


0.00342 



1.8 


0.00249 

0.00306 

0.00173 

1.9 


0.00181 



2.0 


0.00111 

0,00228 

0.00151 

2.2 



0.00172 

0.00157 

1 4 



0.00100 

0.00121 

2.6 



0.00087 

0.00111 

2.8 




0.00085 

3.0 




0.00075 

( e ) 

-r-T/ 2 





TABLE 1 (Continued) 


x/d 



1 

6 

10 

20 

0 

0.00000 

0.00000 

0.00084 

0.00227 

0.1 

0.00000 

0.00000 


0.00132 

0.2 

0.00000 

0.00000 

0.00234 

0.00219 

0.3 

0.00000 

0.00104 



0.4 

0.00000 

0.00130 

0.00207 

0.00520 

0.5 

0.00000 

0.00165 



0.6 

0.00000 

0.00173 

0.00248 

0.00743 

0.7 

0.00014 

0.00177 



0.8 

0.00073 

0.00186 

0.00303 

0.00162 

0.9 

0.00056 

0.00164 

0.00230 


1.0 

0.00036 

0.00129 

0.00210 

0.00266 

1.1 

0.00022 

0.00107 

0.00260 


1.2 

0.00010 

0.00122 

0.00187 

0.00475 

1.3 

0.00005 

0.00141 



1.4 

0.00001 

0.00069 

0.00151 

0.00132 

1.5 

0.00000 

0.00053 



1.6 


0.00463 

0.00138 

0.00213 

1.7 


0.00134 



1.8 


0.00036 

0.00109 

0.00271 

1.9 


0.00033 



2.0 


0.00019 

0.00098 

0.01177 

2.2 



0.00087 

0.00136 

2.4 



0.00060 

0.00139 

2.6 



0.00137 

0.00070 

2.8 




0.00095 

3.0 




0.00131 
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ORIGINAL PAGE sS 
OF POOR QUALITY 


TABLE II 

CONFINED JET FLOWFIELD DATA 


X/D 


R/D 

0.0 

0.5 

1 .O 

1 .5 

2.0 

2 . 5 

0.47500 

0. 13000 

0. 10400 

0. 15300 

0, 15800 

0. 17100 

O . 17900 

0.45000 

0. 14900 

0.11100 

0 , 1 5600 

0. 16500 

0 . 20000 

0 18600 

0.42500 

0. 15800 

0. 1 1300 

0. 16700 

0. 17700 

0.21700 

0 . 20800 

0.40000 

0. 17700 

0 . 1 1 G 00 

0. 18900 

O . 19300 

0.24100 

0 22100 

0.37500 

0. 17700 

0.1 1800 

0.22000 

0.21 600 

0.26500 

0.25000 

0.35000 

0. 18000 

O . 12900 

0 . 26400 

0.25400 

0.29400 

0 27000 

0, 32500 

0. 18500 

0. 16100 

0.32500 

0.28500 

0.34500 

0 . 29400 

0 . 30000 

0. 18800 

0.23000 

0 . 40400 

0.34900 

0.37000 

0. 32500 

0.27500 

0. 19300 

0 . 39000 

0.52000 

0.41800 

0.43300 

O . 35600 

O . 25000 

1 .07000 

O . 63000 

0.65100 

0.47500 

0,49800 

0 39900 

0.22500 

1.02100 

0 . 88500 

0.76400 

0.55 100 

0.53700 

0 42700 

0 . 20000 

1 .00800 

0.97700 

0.85500 

0 . 62 G 00 

0.59800 

0.44900 

0. 17500 

1.00800 

0.99400 

0 . 91 200 

0 . 70600 

0 65900 

0 48 100 

0. 15000 

1 . 00200 

0.98900 

0 . 95000 

0.77300 

0 . 70600 

0.51400 

0. 12500 

0 99900 

0 . 98900 

0 . 96400 

0. 82900 

0.75500 

0.54800 

0. 10000 

1 .00000 

0.98900 

0 . 97000 

0.86200 

0.78800 

0.56300 

0.07500 

0.99700 

0.99100 

0.97200 

0,87600 

0.82400 

0.58500 

0 . 05000 

0.99700 

0 . 99300 

0.98100 

0.88800 

0.84200 

0.59600 

0.02500 

1 .00000 

0.99500 

0.98100 

0.89400 

0 . 84900 

0.62200 

0 . 00000 

1 .00000 

1 . 00000 

0.97800 

0 . 90000 

0 . 86000 

0.63200 


(a) "u/u^ 


ORIGINAL FAr.?- E3 
OF POOR QUALITY 


TABLE II (Continued) 


B/D 

0.0 

0.5 

0.47500 

0.01200 

0 . ooooo 

0.45000 

0 . 00800 

0 . 00300 

0.42300 

0.00300 

0 . 00200 

0.40000 

0 . 00500 

0 , 00500 

0.37500 

0.00300 

0 . 00900 

0.35000 

0.00700 

0 . 00600 

0.32500 

0.00500 

0 . 00600 

0 . 30000 

0.00200 

0 . 0080 C 

0.27500 

0 . 00600 

0.01900 

0.25000 

0.00300 

0.01700 

0.22500 

0.00800 

0.01100 

0 . 20000 

0.00500 

0.00100 

0. 17500 

0.01400 

0 01200 

0. 15000 

0 . 02000 

0.02100 

0 . 12500 

0.02600 

0.02100 

O . 10000 

0.03000 

0 . 02900 

0.07500 

0.02800 

0.03500 

0 . 05000 

0.03300 

0.03700 

0.02500 

0,03000 

O '. 04400 

0 . 00000 

0.03000 

0.04400 


X/D 


1.0 

1.5 

2.0 

0.01 000 

0 . 00500 

O . 00600 

0 . 00800 

0.00700 

0 . 01 500 

0 . 00800 

0.01100 

0.02100 

0.01800 

0.02000 

0.01600 

0 . 01 300 

0,01400 

0.02500 

0.02100 

0.02200 

'\02100 

0.02 100 

0,02400 

0.02200 

0.01600 

0.02600 

0,02400 

0.02500 

0,03100 

C . 03700 

0.01400 

0,02200 

0.03800 

0.01 000 

0.02 100 

0.02900 

0.00100 

0 . 00600 

0,03400 

0 . 00300 

O . 007 00 

0.02500 

0.01100 

0 . 00500 

0.02300 

0.01600 

0.00200 

0.01700 

0.02800 

0.01300 

0 . 01 300 

O . 03000 

0.01900 

0 . 01 300 

0.03900 

0.01900 

0 . 00300 

0.03300 

0,02100 

0 . 00800 

0 . 03600 

0.02400 

0.00400 


2.5 

0 . 00800 
O . 00800 
O O 1 200 
0.01500 
0.01500 
O 01800 
0.01400 
O O^'JCO 
0.02300 
0.01500 
0.02500 
0.02800 
0 . 01 700 
0 0 I IOC 
0 . 00700 
0 . 00300 
0.00100 
0 00700 
0 , 00200 
0.01000 
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TABLE II (Continued) 


X/D 


R/D 

0.0 

0.5 

1 ,0 

1 ,5 

2 .0 

2.5 

0.47500 

0.02700 

0.02600 

0.04200 

0.04500 

0.04900 

0.04900 

0.45000 

0 . 03000 

0.02700 

0.04400 

0.04700 

0.05400 

0 05100 

0.42500 

0.03300 

0.02900 

0.04700 

0.05100 

0 . 059 C 0 

0 . 05600 

0 . 40000 

0.03700 

0 . 02900 

0.05300 

0,05600 

0.06500 

0 . O 6 O 00 

0.37500 

0.03700 

0 . 03000 

0.06200 

0.06300 

0.07200 

0.06700 

0 . 35000 

0 , 03 R 00 

0 . 03400 

0.07300 

0.07200 

0.07800 

0.07300 

0.32500 

0.03900 

0.04400 

0 . 08800 

0 . 08000 

0.09100 

0 . 07800 

0 . 30000 

0.03800 

0 , OG 300 

0, 10200 

0.09400 

0.09800 

0 . 0 B 40 O 

0.27500 

0 . 03900 

0. 10000 

0. 1 1800 

0. 10800 

0. 10600 

0.09100 

0.25000 

0.01400 

0.13000 

0. 13300 

0. 1 1900 

0. 11900 

0. 10200 

0.22500 

0.01000 

0. 10200 

0. 12500 

0.12700 

O . 12400 

0, 10600 

0 . 20000 

0.01000 

0.05600 

0.09900 

O . 12900 

O . 13400 

0 . 1 10 OO 

0 . 1 7500 

0.01000 

0,03100 

0.08000 

0. 12800 

O . 13700 

O . 1 1500 

0. 15000 

O . 00900 

0.02200 

0.06100 

0.11 200 

0. 13400 

0. 12100 

O . 1 2500 

0.01000 

0.01700 

0.05000 

0 . 1 O 10 O 

0 . 1 3000 

0. 12700 

0. 10000 

0 . 00900 

0,01400 

0 , 04000 

0.08700 

0. 12100 

0. 12600 

0.07500 

0.01000 

0.01400 

0.03700 

0.07400 

0 , 11 500 

0. 12500 

0 . 05000 

0.00900 

0.01200 

0 03300 

0.06900 

O . 10500 

O 12400 

0.02500 

0.01000 

0.01200 

0.03300 

0 .06600 

0 . 09900 

0 . 1 2900 

0 . 00000 

0.01000 

0.01 200 

0.03100 

0. 06600 

0 . 09400 

. 0 . 12600 


(c) U» ^/u 
rme o 
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TABLE II (Continued) 


x/o 


R/D 

0.0 

O.S 

1 0 

1 .5 

2.0 

2.5 

0.47500 

0.01200 

0.01300 

0. 02 800 

0.02800 

0.02800 

0.03000 

0.45000 

0.01300 

0.01600 

0.03100 

0.03200 

0.03600 

O 03600 

0.42500 

0.01500 

0.01700 

O . C 3500 

0.03600 

0.04200 

0.04300 

0 . 40000 

0.01600 

0. 01800 

0. 04000 

0 , 04 100 

0,04900 

0.04800 

0.37500 

0.01700 

0.01900 

0.04800 

0.04700 

0.05500 

0 , 05400 

0 . 35000 

0.01700 

0.02300 

■' .05800 

0.05600 

0 . 06200 

0 06000 

0 . 32500 

0.01800 

0.03200 

0.07000 

0. 06200 

0 , 072.00 

0.06400 

0 . 30000 

0.01700 

0 . 05000 

0.08300 

0.07400 

0.07900 

0.07200 

0.27500 

0.01600 

0.08200 

0.09500 

0.08400 

0.08900 

0.07800 

0.25000 

0 . 00800 

0. 10000 

0.09700 

0 . 09000 

0.09700 

0 . 08400 

0.22500 

0.00600 

0 . 07900 

0 . 09000 

0. 09 300 

0. 10100 

0 . 08700 

0 . 20000 

0.00700 

0.04300 

0 . 07 K . 3 O 

0.00920 

0. 10300 

0.09000 

0. 17500 

0.00600 

0.02400 

0.05400 

0.08800 

0. 10400 

0 . 09400 

0. 15000 

0.00700 

0.01500 

0.04200 

0 . 08000 

0. 10200 

0 . 09500 

0 . 12500 

0.00600 

0.01 100 

0.03100 

0.06900 

0.09600 

O . 09500 

0 . 10000 

0.00600 

0.00900 

0.02500 

0.06000 

0.08800 

0 . 09600 

0.07500 

0.00600 

0 . 00700 

0,02200 

0.05300 

0.08300 

0.09500 

0.05000 

0.00600 

0.00700 

0.01900 

0.04600 

0. 07800 

0.09300 

0.02500 

0.00600 

0 . 00700 

0.01600 

0.04400 

0.07300 

0 . 09600 

0 . 00000 

0.00700 

0.00700 

0.01800 

0.04400 

0.07200 

0.09700 


/u 

nns o 


( d ) V * 
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TABLE II (Continued) 


X/D 


R/D 

0.0 

0.5 

1 .o 

1 5 

2.0 

2 . 5 

0.47500 

0.00000 

0 . 00000 

0.00029 

0.00025 

0.00033 

0.00045 

0 . 45000 

0 , 00000 

0 . 00002 

0 . 00033 

0.00032 

0.00054 

0.00049 

0.42500 

0.00002 

0 . 00004 

0 . 00043 

0 . 00038 

0 . 00070 

0.00072 

0.40000 

0 . 00000 

0.00008 

0 00063 

0 . 00055 

0.00104 

0.00093 

0.37500 

0.00001 

0.00010 

0.00100 

0.00084 

0.00142 

0.00136 

0 . 35000 

0 . 00003 

0 . 0002 1 

0.00161 

0.00126 

0.00197 

0.00172 

0. 32500 

0.00007 

O . CK )063 

0.00278 

0.00183 

0.00284 

0.00216 

0.30000 

0.00010 

0.00154 

0 . 004 B 6 

0 . 003 1 2 

0.00376 

0 . 0028 1 

0.27500 

0.00117 

0.00513 

0.00839 

0.00466 

0.00521 

0.00335 

0.25000 

0.001 12 

0.01252 

0.01272 

0.00645 

0.00694 

O 00431 

0. 22500 

0.00029 

0.01703 

0.01582 

0.00879 

0.00850 

0.00484 

0 . 20000 

0 . 000 1 3 

0.00239 

0.01443 

0.01 134 

0.01038 

0 00527 

0 . 17500 

0 . 000 1 1 

0 . 00004 

O 01021 

0 . 01 308 

0.013 to 

0.00557 

0. 15000 

0.00007 

0 . 00004 

0.00539 

0.01353 

0.01330 

0 . 00595 

0. 12500 

0 . 00007 

0.00006 

0.00056 

0.01226 

0.01410 

0.00631 

0. 10000 

0 . 00006 

0 . 00007 

0.00032 

0 . 00948 

0.01367 

0.00630 

0.07500 

0.00006 

0 . 00008 

0 . 00036 

0.00769 

0.01255 

0 00547 

0 . 05000 

0.00005 

0 . 00006 

0 . 0004 1 

0 . 0055,1 

0.01121 

0.00459 

0.02500 

0.00002 

0 . 00007 

0 . 00054 

0.00499 

0.00942 

0.00374 

0 . 00000 

0.00002 

0.00011 

0.00073 

0.00487 

0.00792 

0.00202 


(e) u*v’/u^^ 

o 
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TABLE II (Continued) 


R/D 

0.0 

0.5 

X/D 

1 ,o 

1 .5 

2.0 

2.5 

0,47500 

0 . 00000 

0 . OOfJ 0 1 

0. 00490 

0.00182 

0 . 00058 

0.00326 

0 . 4500 C 


0.00042 

0,00239 

0.00169 

0.001 10 

0.00172 

0,42500 

0.00001 

0 . 00083 

0.00131 

0.00135 

0.00172 

0 . 00207 

0,40000 

0.00003 

0.00157 

0,00120 

0.00143 

0.00219 

0.00224 

0,37500 

0.00034 

0 . 00074 

0.00134 

0.00138 

0.00270 

0.00279 

0,35000 

0 . 00039 

0 . 00049 

0,00155 

0.00184 

0,00248 

0 . r 0395 

0,32500 

0.00003 

0.00063 

0.00200 

0.00194 

0.00377 

0.00397 

0,30000 

0 . 00097 

0 . 00068 

0 . 0025 1 

0.00237 

0.00431 

0 . O 0457 

0.27500 

0.00014 

0.00129 

0.00342 

0.00373 

0.0044 1 

0.00456 

0.25000 

0.00013 

0.00255 

0.00525 

0.00489 

0.00672 

0 . 006 1 2 

0,22500 

0,00047 

0.00495 

0.00782 

0 00587 

0.00786 

0.00976 

0 , 20000 

0.00107 

0 . 0022 1 

0.00983 

0.00738 

0.00857 

0.00983 

0. 17500 

0.00181 

0.00034 

0.01083 

0.00897 

0.01129 

0.00864 

0. 15000 

0.00079 

0 . 0008 1 

0.01045 

0.01 109 

0.01396 

0 00896 

0. 12500 

0.00348 

0.00000 

0.00282 

0.01388 

0.01734 

0.01299 

0, 10000 

O . 00304 

0 . 0035 1 

0.00397 

0.02031 

0.01997 

0.01716 

0,07500 

0.00154 

0.00191 

0.00333 

0.02985 

0.02343 

0.01670 

0.05000 

0 . 00099 

0.00159 

0.00465 

0.03087 

0.04526 

0.01250 

0.02500 

0 . 00000 

0.00107 

0. 01823 

0.04 191 

0.05261 

0.01047 

0.00000 

0.00000 

0 . 00000 

0 . 00000 

0 . 00000 

O . 00030 

0 . 00000 


(f) V^* 
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TABLE III 

CONFINED JET WITH CONTRACTION 




NOZZLE 

FLOWriELD 

DATA 


R /0 

0.0 

0.5 

X/D 
1 .0 

1 .5 

2.0 

0.47500 

0. 15700 

0. 10100 

0. 15400 

0. 18700 

0.22600 

0.45000 

0. 16500 

0. 10400 

0 . 1 6400 

0. 19600 

0.22700 

0.42500 

0. 16600 

0. 10700 

0 . 17100 

0.21400 

0.24100 

0 . 40000 

0. 18400 

0. 1 1000 

0. 18900 

0.23800 

0.25700 

0.37500 

0. 19500 

0 . 1 2000 

0.22800 

0.27200 

0.29300 

0.35000 

0. 19400 

0. 13900 

0.27300 

0 30500 

0. 31400 

0.32500 

0 . 20300 

0. 19000 

0.32800 

0.35100 

0 . 36000 

0 . 30000 

0.2 1400 

0.29100 

0.40300 

0.42500 

0.41 300 

0.27500 

0.21000 

0.45600 

0.51300 

0 . 50400 

0.47400 

0.25000 

0.20700 

0.71400 

0.63100 

0.57600 

0.52300 

0.22500 

1 .02400 

0.91800 

0.75100 

0.64500 

0.58000 

0 . 20000 

0.99600 

0.98700 

0.85000 

0 .73700 

0.61800 

0. 17500 

0 - 99900 

0 99000 

0.92300 

0 . 78800 

0.67600 

0 . 15000 

0.99300 

0.99300 

0.95800 

0.83800 

0.73300 

0 . 12500 

0.98400 

0.99300 

0 . 97800 

0.87200 

0.76600 

0 . 10000 

0.98700 

0.99300 

0 98600 

0.89400 

0.79800 

0.07500 

0.98800 

0 . 99600 

0 . 98900 

0 . 90600 

0 . 80900 

0.05000 

0.99 100 

O . 99400 

0.98900 

0.92000 

0.82500 

0.02500 

0 . 98800 

0 . 996 C 0 

0 .99600 

0.92600 

0.83500 

0.00000 

0.99700 

1.00000 

0.99100 

0 . 92600 

0 . 84000 


(a) u/u^ 


TABLE III (Continued) 


x/o 


R/D 

0.0 

O.S 

1 ,o 

1 ,5 

2.0 

0.47500 

0 . 00200 

0.00900 

0 . 00700 

0 . 00800 

0.00900 

0.45000 

0 . 00200 

0.00700 

0 . 00800 

0 . 00800 

0 . 00200 

0.4 2500 

0 . 00200 

0 . 00800 

0 . 00600 

0.01400 

0.01 1 CK - 

0 . 40000 

0.00400 

0 . 00500 

0.00200 

0.02700 

0.01600 

0.37500 

0.00400 

0 . 00300 

0 , 00400 

0 . 01 500 

0.01800 

0 , 350 GO 

0.00700 

0.00100 

0 . 00500 

0.02200 

0.01700 

0 . 32500 

0 . 00800 

0 . 00000 

0 . 00800 

0 . 03800 

0 . 02000 

0 , 30">00 

0 . 00400 

0 . 00600 

0 . 00800 

0.03100 

0 . 01 400 

0.27500 

0 . 00200 

0.00700 

0.01500 

0.02900 

0,01200 

0 . 25000 

0 . 00800 

0 00500 

0.01700 

O . 03600 

O . 00800 

0.22500 

0.01 200 

0.00700 

O . 00300 

0 . 02900 

0 . 00500 

0 . 20000 

O . 00600 

0.01500 

0.01 100 

0.02400 

0 . 00900 

0 , 17500 

0,02100 

0 . 02400 

0.01000 

0.00700 

0 . 00300 

0 . 15000 

0.02700 

0 . 03000 

0 , 02000 

0.01000 

0.01 100 

0 . 12500 

0.03600 

0 . 0340 C 

0.02700 

0 . 00000 

0 . 01 300 

0. 10000 

0.04300 

0 . 039 CI 0 

0.03200 

0 . 00900 

0 , 00900 

O . 07500 

0.04500 

0.04500 

0 . C 3800 

O . 00400 

0 . 00200 

0 . 05000 

0.04800 

0.04300 

0.04700 

0.01800 

0 . 01 400 

0.02500 

0.04500 

0.04500 

0.04500 

0.02500 

0.01200 

0.00000 

0.04600 

0 . 04600 

0.05400 

0 . 02500 

0.02200 


(b) v/u 
o 
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TABLE III (Continued) 


R/D 

0.0 

0.5 

X/D 
1 .O 

1 .5 

2.0 

0.47500 

0 . 02900 

0.02300 

0.03900 

0.05200 

0.06500 

0.45000 

0.03100 

0.02300 

0.04300 

0.05400 

0.06700 

0.42500 

0.03100 

0.02400 

0.04400 

O . 06000 

0.06900 

0.40000 

0.03600 

0.02500 

O . 05100 

O 06700 

0.00730 

0.37500 

0.03800 

0.02900 

0 . 06000 

0.07500 

0.08300 

0.35000 

0.03800 

O . 03 SO 0 

0.07100 

0 . 082 C 0 

0.08700 

O . 32 SO 0 

0.03900 

O . 05000 

0.08300 

0 . 09500 

0,09700 

0 . 30000 

0 . 04 100 

0.07800 

O . 09900 

0. 1 1000 

0. 10900 

0.27500 

0.04100 

0. 11500 

0 . 1 1900 

0. 12300 

0 . 1 2000 

0.25000 

0,03100 

0. 14000 

0 . 1 3600 

0. 13800 

O . 12700 

0.22500 

0.01300 

0. 10000 

0. 13300 

0. 13800 

0. 13400 

0 . 20000 

0 . 0)200 

0 . 05500 

0 , 12)00 

0. 14400 

0. 13600 

0 . 17500 

0.01200 

0 . 03500 

0.09500 

0. 12600 

0. 13500 

0 , 15000 

0.01300 

0.02500 

O , 07600 

0 . t 1400 

0. 13400 

0 . 12500 

0.01200 

0 . 02000 

0 . 06000 

0 09900 

0, 12600 

0. 10000 

0.01300 

0.01800 

0,05100 

0.08400 

O . 12300 

0.07500 

0.01200 

0.01700 

O . 04400 

0.07900 

O . 1 1300 

0.05000 

0.01300 

0.01600 

0.04300 

0.07000 

0. 1 1200 

0.02500 

0.01200 

0.01500 

0 . 03900 

0 . 06600 

0 10700 

0.00000 

0.01300 

0.01600 

0.03700 

0 06400 

0. '0800 


(c) u* /u 
' ' nns 0 
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TABLE III (Contlnucid) 


x/o 


R/p 

0.0 

0.5 

1.0 

1 . 5 

2.0 

0.47500 

0.01200 

0.01200 

0.02500 

0.03100 

0.03600 

0.45000 

0.01300 

0. 01400 

0 . 03000 

0-03700 

0.04000 

0.43500 

0 . 01 300 

0.01500 

0.03300 

0.04300 

0.04500 

0 . 40000 

0.01500 

0.0 1600 

0.03800 

0 . 05000 

0 . 05000 

0 . 37500 

0.01500 

0.01800 

0.04600 

0.05800 

0.05900 

0.31000 

0.01500 

0.02400 

0.05800 

0.06500 

0 . 06600 

0 . 32500 

0.01600 

0.03700 

0.07100 

0.07400 

0.07400 

0 . 30000 

0 . 0 I 600 

0.06400 

0 . 08600 

0 08700 

0.08500 

0.27500 

0. 01400 

0.09900 

0. 10300 

0 09700 

0.09400 

0.25000 

0 . 00800 

0. 11500 

0. 1 1300 

0. 10100 

0. 10000 

0.22500 

0 . 00700 

0.08200 

0. 10900 

0. 10000 

0. 10400 

0,20000 

0 . 00600 

0.04500 

0.09300 

0 . 09800 

0. 10300 

0 . 17500 

0.00600 

0.02500 

0.07100 

0.08500 

0. 10100 

O . 15000 

O . 00600 

0. 01700 

0.05300 

0.07600 

0.09700 

0. 12500 

0.00600 

0.01200 

0.03900 

0.06100 

0.09400 

0.10000 

0.00600 

0.00900 

0.03200 

0 . 05600 

0 08600 

0.07500 

0.00600 

0.00800 

0 . 02 G 00 

0.04800 

0.07900 

0 . 05000 

0.00600 

0.00800 

0.02300 

0.04200 

0.07700 

0.02500 

0 . 00600 

0 . 00700 

0 . 02000 

0.04100 

0.07600 

0.00000 

0. 00600 

0 . 00700 

0.01900 

0.03900 

0.07500 


(d) V* /u 

' ■' nos o 
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TABLE III (Continued) 


R/D 

0.0 

0.5 

X/D 
1 .0 

1,5 

2.0 

0.47500 

0.00000 

0 . 00008 

0.00019 

0.00029 

0 . 00039 

0 . 45000 

0 . 00000 

0.00009 

0 . 00033 

0.00045 

0.00046 

0.42500 

0.00001 

0.00015 

0,00043 

0 . 00068 

0.00069 

0 . 40000 

0.00001 

0 . 00020 

0.00062 

0.001 10 

0.00095 

0.37500 

0.00001 

0 , 00026 

0.00102 

0.00159 

0.00164 

0.35000 

0.00001 

0.00039 

0.00169 

0.00240 

O 00226 

0 . 32500 

0.00003 

0.00090 

0.00272 

0.00337 

0.00350 

0.30000 

0 . 00003 

0.00271 

0.00449 

0.00548 

0 . 005 1 1 

0.27500 

0.00007 

0 , 00707 

0,00763 

0 . 008 1 2 

0.00731 

0.25000 

0.00312 

0,01594 

0.01 132 

0.01056 

0.00927 

0.22500 

0 . 00006 

0.01476 

0.01422 

0,01287 

0.01 100 

0 . 20000 

0,00005 

0.00026 

0,01431 

0.01507 

0.01270 

0. 17500 

0 . 00005 

0 , 00004 

0.00150 

0.01471 

0.01410 

0. 15000 

0.00005 

0.00004 

0 . 00052 

0.01332 

0.01419 

0 . 12500 

0.00004 

0.00001 

0 . 00005 

0 . 005 1 1 

0.01407 

0. 10000 

0 . 00004 

0 . 00004 

0 , 00002 

0.00475 

0.01331 

0.07500 

0.00004 

0 . 00006 

0 . 00002 

0.00314 

0.01155 

0.05000 

0.00003 

0 . 00005 

0 . 00005 

0.00266 

0.01035 

0 . 02500 

0.00003 

0.00005 

0 . 00006 

0.00298 

0 . 00900 

0.00000 

0 . 00003 

0 . 00005 

0 , 000 1 1 

0.00297 

0.00720 

(D 

p » 2 
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TABLE III (Continued) 


R/D 

0.0 

0.5 

X /0 
1 .0 

1 ,5 

2.0 

0.47500 

0.00002 

0.00140 

0 . 00098 

0.00164 

0.01878 

0.4 5000 

0 , OOOCF 

0.00145 

0.00195 

0.00167 

0 . 00307 

0.42500 

0.00003 

0.00254 

0.00175 

0.00163 

0.00233 

0 . 40000 

0 . 00003 

0.00151 

0 .O 0109 

0.00191 

0.00136 

0.37500 

0 . 00007 

0 . 0009 1 

0.00123 

0.00240 

0.00223 

0.35000 

0,00009 

0.00057 

0.00170 

0 . 00307 

0.00331 

0.32500 

0,00015 

0.00050 

0 . 002 1 1 

0.00283 

0.00368 

0 . 30000 

0 . 00042 

0.00103 

0.00245 

O . 0036 1 

0.00486 

0.27500 

0 . 00095 

0.00168 

0.00337 

0.00542 

0 , 007 1 5 

0.25000 

0,00039 

0.00348 

0.00480 

0 00756 

0.00935 

0.22500 

0.00001 

0.00545 

0.00654 

0 , 00806 

0.01260 

0 . 20000 

0 . 0002 1 

0,00036 

0.00839 

0.01063 

0.01376 

0 . 1 7500 

0.00182 

0 . 00070 

0.00140 

0.01468 

0.01367 

0 , 15000 

0.00033 

0.00126 

0 . 00095 

0.01598 

0.01932 

0. 12500 

0.00074 

0 . 00000 

0.00016 

0 . 00920 

0.02086 

0 . 10000 

0.00109 

0.00149 

0.00020 

0.01407 

0.03629 

0.07500 

0 . 0002 S 

0 . 00594 

0 . 00058 

0.01219 

0.05544 

0 . OSOOO 

0 . 00000 

0 . 00000 

0.00075 

0,01347 

0.03157 

0.02500 

0 . 00053 

0.00078 

0.00284 

0.05014 

0 . 00000 

0.00000 

0,00018 

0.00065 

0.00109 

0 . 00000 

0.00070 


<f) V^* 
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Figure 3, The Confined Jet Flowfield 





Figure 4. Free Jet Apparatus 
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LM741CN Operational Amplifier 


Adder Circuit 


100 KJX 



Subtractor Circuit 

Figure 9. Schematic of the Constructed Devices 
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Figure 10. Instrumentation Layout 
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Figure 14, Angulation Test Coordinate System 





















SUIJ 










APPHNDIX C 


COMPUTER PROCRAM USER'S GUIDE 


sre?!^v- 


85 


A listing for the computer program discussed in this Appendix and 
sample addressed data files are contained In Appendix Included it 
the end of this user's guide are the nomenclature employed in the pro- 
gram and a sample output. 

The X“Wire arrangement used for data input into this code is assum- 
ed to be the standard mutually perpendicular, 45° to probe axis wire 
orientation. The hot-wires are assumed to obey an extension of King's 
Law. That is: 

2 — ^ — 

E = A + Bu + Cu 

where. A, B, and C are constants determined a prio ri by a calibration 
experiment. 

Four main data sets are needed as input to generate the data reduc- 
tion. They are: 

1. Calibrnt'on data 

2. Expertmer;^:al dat 

3. Miscellaneous data consisting of the date, experiment 
parameters and ambient conditions. 

4. Angulation data 
and are now itemized in detail. 

Calibration 


The calibration is accomplished as discussed in section 3.1. Thus, 
the calibration data set contains a number (NMAX) of selected velocities 
and the corresponding mean voltages from wire 1, CALEl and wire 2, CALE2. 
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Experimental 

This data 1 b the Information from the test section. At each radial 
position, r/d or r/R, seven quantities are measured. They are: 

1. Mean voltage of wire 1, 

2. Mean voltage of wire 2, 

3. Adder voltage (ej^' + ADDE 

4. Subtractor voltage (e^^* ~ ®2*^’ SUBE 

5. Attenuation of channel A of the correlator, ATTA 

6. Attenuation of channel B of the correlator, ATTB 

7. Correlator output voltage (millivolts, DC), EMVDC 
There are jmax of each seven quantities. 

Mlacellaneous 

The mlscellantous data set is just that, miscellaneous. It con- 
tains : 

1. Month, day, year - the date of the experiment 

2. XDD - the downstream location of the probe in diameters 

3. LDD “ the downstream position of the contraction block in 
diameters 

4. ALPHA - the expansion block angle in degrees 

5. FE - the swirl angle, should equal zero 

6. PAMB - the ambient pressure, cm of hg 

7. TAMB - the ambient temperature, “C 

8. E«>1 - nozzle exit velocity voltage of wire 1, volts 

9. E02 - nozzle exit velocity voltag j of wire 2, volts 

10. EOl ~ no flow voltage of wire 1, volts 

11. E02 - no flow voltage of wire 2, volts 




C 
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Angulation 

Because the angulation test Is only done once» the slopes of and 
E 2 versus yaw angle, <j) are not calculated inside this computer code. 
Iherefore, this data consists of a series of the two slopes 3Ej^/ 34> and 
and their corresponding velocities. 

Computer Code 

This computer code consists of a main program and three subroutines. 
Bach is separated into chapters and sections. Thus, itemized by 
chapters : 

Main 

1. Inputing Data 

2. Curve Fitting the Calibration Data (calls CFIT) 

3. Calculating the Sensitivity Coefficients 

4. Curve Fitting D and B (calls C2FIT) 

m 

5. Calculating the Jet Nozzle Exit Velocity 

6. Calculating All Dutput Variables 

7. Data Output in Tabular Form 

8. Data Output in Graphic Form (calls SUBPLT) 

Subroutines CFIT and C2FIT 

1. Initializing the Variables 

2. Repackaging the Curv€>. Fitting Variables 

3. Calculation of A, B, and C, Calibration Constants, 

In CFIT and Slope, B and Y-intercept, A in C2FIT 
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SubroutlnR Subpit 
1. Line i’lot Graph Routine 

The following Is an explanation of each chapter and section. Refer to 
Appendix D for a code listing. Nomenclature for the computer code Is at 
the end of this appendix. 

Main 

Chapter 1 

Section 1 

This section Inputs the code parameters of the number of velocities 
used in calibration; NMAX, the number of velocities used la the angula- 
tion tests; IMAX* and the number of measuring positions used on a cross 
section In the test section experiment, JMAX. 

Section 2 

CALU is Inputed through a data statement because the calibration of 
the probe is done at these velocities every time. The data of the cali- 
bration voltages from the hot-wires Is read from a separate data file 
which is called Into the code using statement 4240. 

Sectltm 3 

Like GALU, RDD doesn't vary. Thus, RDI) Is contained in a data 
statement. The seven quantities measured at each location in the test 
section are entered in a data set. As before, this data Is called in 
the code using statement 4250. 
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Section 4 

The data from the angulation tests Is contained In this section. 
Data statements are used as the Input medium because these tests are 
done once per probe. 

Section 5 

Miscellaneous data is Input hetre through a allocation statement 
that addresses the miscellaneous data file (statement 4260). 

Chapter 2 

Section 1 

The calibration data for wire 1 from Chapter 1, Sec:tlon 2, is curve 
fit using a second order least squares routine (CFIT) sc that: 

2 3 - 

CALEl = A1 + Blu + Clu 


Section 2 


The calibration data for wire 2 Is curve fit as wire 1 so that; 


CALE2^ 


H 

A2 + B2u + C2u 


Chapter 3 

In this chapter^ I> and B, are calculated for the velocities given 
m 9 

In the angulation experiment. 
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Chapter 4 


Section 1 


The D 's of the previous chapter are linearly curve fit so that: 


D = A3 + B2 u 
m 


Section 2 


The B.*s of the prevloxis chapter .*re linearly curve fit so that: 

B . o A4 + B4 u 
<P 

Chapter 5 

The nozzle exit velocity, , used to nondlmenslonallze all turbu- 
lent quantities is calculated here from two hot-wire mean voltages 
determined from positioning the probe at the uozzle exit face. 


Chapter 6 

u'/u , v'/u and u'^v’/u ^ are calculated, 
o o o 

Chapter 7 

Section 1 

Miscellaneous data is written out here with appropriate headings. 
Section 2 

This section writes the headings for the major output. 

Section 3 

Section 3 writes the turbulent quantities In tabtilar form. 


The outputs; ^/u , ^/u , 
o o 
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Chapter 8 

Hils chapter then graphs each turbulent quantity using a line printer 
plotting routine that Is a library function in the machine for which this 
code was written. To do the plots, a subroutine (SUBPLT) Is called to 
access the library plotting function. 

Subroutine CFIT 

Chapter 1 

The variables are Initialized to zeio. 


Chapter 2 


The X and y are repackaged so that the subroutine curve fits CALEl 

2 

versus u the first time It is called and CALE2 versus u the second. 


Chapter J 

The necessary calculation are done using the least square technique. 
Ihe last part of this chapter uses Cramer's Rule to solve for the A, B, 
and C constants. 

Subroutine C2FIT 

Chapter 1 

The variables are Initialized to zero. 


Chapter 2 

The X and y are repackaged so that the subroutine curve fits D 

versus u the first time it Is called and B. versus u the second. 

<P 




s 


Chapter 3 


A linear curve fit la applied to D and veraua u to calculate 

in 0 

slopes and y-lntercepta. 


Implementation 

To use this computer code, the following cards have to be changed. 

1. NMAX (line 380) 

2. IMAX (line 390) 

3. JMAX (line 400) 

4. Line 4240 will be changed to address the file of calibra- 
tion data that corresponds to the particular experiment. 

5. Line 4250 will be changed to address the file of experi- 
ment data that corresponds with the particular experiment 

6. Line 4260 will be changed to address the file of miscel- 
laneous data that corresponds to the particular experi- 
ment. 

7. If one has crossed-wlres that produce different 3E/3^'s, 
then lines 800 and 810 will have to be modified. 

It must be remembered that the code has the following restrictions: 

1, Non-swirling flow 

2, Wires obey King's Law 

3, Wires are matched so that (e^^' + e^*) and (e^^' - e^*) 
can be simplified. 


4, The hot-wire geometry Is the same as shown In Figure 7 
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NOMENCLATURE 


Main 


Code 

llieory 

Meaning 

ADDE 

(ej- + .j’) 

Adder voltage, volts 

ALPHA 

a 

Wall expension angle, degrees 

AMI 

A 

“l 

Axial velocity fluctuation sensitivity 


coefficient of wire 1 

AM2 

A 

Axial velocity fluctuation sensitivity 



coefficient of wire 2 

ANGU 


Angulation experiment velocity, iii/sec 

ATTA 

ATTA 

Attenuation of channel A 

ATTB 

ATTB 

Attenuation of channel B 

AVI 

A 

Vj 

Radial velocity fluctuation sensitivity 


X 

coefficient of wire 1 

AV2 

A 

Radial velocity fluctuation sensitivity 


1 , 

coefficient of wire 2 

Al, A2, etc. 

^ 2 * etc. 

Y~intercept of linear curve fit 

BO 

®<t> 

(E A + E, A ) 

1 Vl 2 Vj 

Bl, B2, etc. 

®2* 

Slope of linear curve fit 

CALEl 

Calel 

Calibration voltage of wire 1, volts 

CALE2 

Gale 2 

Calibration voltage of wire 2, volts 

CALU 

- 

Calibration velocity, m/sec 

DM 

D 

tn 

(E, A + E. A ) 
1 ^ ®2 

TAMB 

TAMB 

Temperature of atmosphere, °C 

U 

- 

Axial velocity, m/sec 


Average of nondinenslonal axial velocity 


UAVG 
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UPDUO 

u'/u 

o 

Nondlmenslonallzed rms axial velocity 



fluctuation 

UPRMS 

u* 

SMS axial velocity fluctuation, u'rms, 



m/sec 

UO 

u 

0 

Nozzle exit velocity, m/sec 

U1 


Time mean axial velocity of wire 1, m/ '-^c 

U2 

“2 

Tlnffi mecii axial velocity of wire 2, m/sec 

UIDUO 

c 

0 

Nondlmensional axial velocity, wire 1 

U2DUO 

“2'“o 

Nondlmeneional axial velocity, wire 2 

VPDUO 

v'/u 

0 

Nondlmenslonallzed rms radial velocity 



fluctuation 

VPRIffi 

V» 

RMS radial velocity fluctuation, v’rms, 



m/sec 

XDD 

x/d or x/D 

Nondlmenslonallzed distance from jet exit 

XDL 

x/L 

Nondlmenslonallzed position of constriction 



block 

E 

“ 

Variable defined by call statement 

EMVDC 

EMVDC 

Correlator output, millivolts, D.C. 

EOl 

- 

Voltage of wire 1; no flow 

E02 

- 

Voltage of wire 2; no flow 

El 


^tean voltage of wire 1, volts 

E2 


Mean voltage of wire 2, volts 

EIDFE 

3E^/3(j) 

The slope of E^^ versus (|i, volts/rad 

E2DFE 

3E2/3«^ 

The slope of E^ versus <|i, volts/rad 

E1D2U 

3Ej^/3u 

The slope of versus u, volts/m/sec 

E2D2U 

3 E 2 / 3 U 

The slope of E 2 versus u, volts/m/sec 
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EOl 


Nozzle exit velocity voltage of wire 1, 
VOltfi 

E02 


Nozzle exit velocity voltage of wire 2, 
volts 

FE 


Swirl angle (zero for this program), degrees 

IMAX 

IMAX 

Number of points in angulation experiment 

JMAX 

JMAX 

Number of points in test section experiment 

NMAX 

NMAX 

Number of points in calibration experiment 

PAMB 

PAMB 

Pressure of atmosphere, psla 

REYSTR 

u*v' 

=~r 

Nondlmenslonal rx component of Reynolds 


% 

stress tensor 

RDD 

r/d or r/D 

Nondlraensionallzed distance from center line 

SUBE 

Cl* - e^') 

Sub tractor voltage, volts 

Subroutine CFIT 


AAO, AAl, etc* 

•• 

Constants for Cramer’s rule 

A, B, etc. 

- 

Curve fit variables 

BI, Cl, etc. 

-- 

Determinate constants 

DET 

- 

Value of determinate used In Cramer’s rule 

E 

- 

Voltage 

SI 

- 

Value of NMAX 

SX 

- 

Sum of X 

SXY 

- 

Sum of X times y 

SX2 

- 

Sum of u 

SX3 

- 

Sum of u times x 

SX4 

- 

2 

Sum of u 

SX2Y 

_ 

Sum of x^ times y 
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SY 

U 

X 

X2 

X3 

X4 

XY 

X2Y 

Y 


2 

Sum of E 

Time-mean axial velocity 

li 

Abclssa value - u 
Dummy variable 

u times x 

2 

u 

X times y 
2 , 

X times y 

7 

Ordinate value - E 


Subroutine C2FIT 

ANGU 

CONST 

SX 

SXY 

SX2 

SY 

U 

X 

XMEAN 

Y 

YMEAN 


Angulation test velocities 

Variable defined by call statement 

Sum of X 

Sum of X times y 
2 

Sum of X 
Sum of y 

Variable defined by call statement 

Abclssa axis values 

Mean value of abclssa values 

Ordinate axis values 

Mean value of ordinate values 


Subroutfue SUBPLOT 

IM - Field size 

NLAST - Number of points to be plotted 

Abclssa plot variable 


X 
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OR1QO?JAE. n 

OF POOR QUALilV 


SAMPLE OUTPUT 


OAtA TAKtN »/a3/e3 



K/Oi 1 . ft 
AiPHAtOe . 0 
POM0I74.3 TAMDi 

1 27 0 

U0< 10 3430 

l/P* 2 0 
ro • 0.0 
eoi *2 01 

e02i 2 .07 


fl/0 

Ut/UO 

U2/UO 

U/UO 

v/uo 

U“/UO 

V- /uo 

U“ V-/U0*»2 

0 . 000 

O . Oft 1 

o.oet 

0 024 

O 034 

0 004 

0 034 

0 002730 

0 02S 

0 Oft 1 

0.001 

O 024 

0 03ft 

0 ooo 

0 034 

0 00374 1 

0 . o*r ; 

0 »34 

0 OO I 

O 020 

0 014 

0 070 

0 030 

0.002440 

. 

0 . ft lo 

0.001 

Q . 004 

e 004 

0 070 

0 046 

0.002000 

^ fOO 

0 003 

0 . 000 

O 404 

0 OOfl 

0 004 

0 Oftt 

0.004302 

0 120 

0 ft72 

0 072 

O 072 

0 . ooo 

0 OOO 

0 OftO 

0 . 004 607 

0 ISO 

e . 024 

O .044 

0 034 

0 010 

O 14 4 

0.070 

0 012234 

0.17ft 

0 7ft 1 

0.700 

0 704 

0 007 

O 136 

0 OTO 

0 . 013ft 1 1 

0 . 200 

0.713 

0.701 

0 737 

0 034 

0.144 

0 ooo 

0 013440 

0 22ft 

0 oto 

O . B?3 

0 04ft 

0 02ft 

0 134 

0 002 

0 . 0114 1 1 

0 26^ 

0.041 

0 012 

0 074 

0 036 

0 1 SO 

O 002 

0 . 000700 

0 27ft 

0 4?: 

0 . B33 

0 ft04 

0 02ft 

0 1 23 

0 OOft 

0 . 007446 

C . 300 

0 . 304 

0 4B7 

0 .424 

0 03 1 

0 no 

0 000 

0 . 000033 

0 S2ft 

0.314 

0.300 

0 301 

0.034 

O Oftft 

0 . 004 

0 . 003000 

0. 3ftO 

O . 202 

O. 327 

0 306 

0 022 

0.007 

0 ooo 

0.002202 

0 . 37ft 

O . 3B8 

0 .307 

0 272 

0.014 

0 07ft 

0 0B3 

0.001440 

0 . 400 

0 . 2 M 

0 .204 

0 . 2ib 

0 . 027 

0.007 

0 . 046 

0 . 00100B 

0 . 42ft 

0.201 

0 . 220 

0.214 

0.014 

0.000 

e 030 

0.000633 

0 . 4ftO 

O too 

0 . 204 

0 114 

O 004 

0 Oft« 

O 034 

0 . 000410 

0 . 47S 

O. 170 

0 . 10ft 

0 147 

0 OOft 

0 . 0B2 

0 . 020 

0 . 000300 





n 




1 

1 . 


I 


l 

il 

i 


t 

f 


{ 




ORIGINAL PAO!£ U 
OF POOR QUAUIY 

■ /o 

u/uo 

^ 0^ 0 

o 634 

O.Oflft 

0 .634 

0 060 

0 .630 

0,098 

0 . 606 

0 . 100 

0-664 

0. 13« 

O 693 

0 . lie 

0 . 096 

0 . 196 

0 960 

0-300 

0 .939 

0 89B 

n . 646 

O.360 

0,694 

0 998 

0 . 004 

o . noo 

o 436 

0.936 

0.961 

0 . 360 

0 . 906 

0 . 37B 

0.393 

0 . 400 

e . 334 

0 436 

0.314 

O. 460 

0. 144 

O 496 

0 169 


0.000 

0.000 


• 600 
U/l/0 





1 

■j 


\ 


i 



ORIGINAL PAGE IS 
OF POOR QUALITY 


A/D 

7/00 

0 000 

0 02B 

0 029 

0.029 

O 090 

O .Ot« 

0 079 

0.004 

0 too 

0 . 009 

0 129 

0 . ooo 

O 190 

0.010 

O. 179 

0,007 

0 . 200 

O . 024 

0 . 229 

O . 029 

O 2«0 

0.039 

0.279 

0.029 

0 . 300 

0.03 1 

0 . 326 

0 031 

0 . 390 

0 . 029 

0 . 379 

0.016 

O . AOD 

0.027 

O 939 

0.014 

0 I90 

0 . 009 

0 479 

0 . 009 


i 









ORlGWi^'^- ^ 

OF POOR QUftLOT 


«/0 

U*/U0 

O . 000 

0.004 

o. oaft 

0 .000 

0.000 

0.070 

0 OYl 

0.070 

0.100 

0.004 

0. ia« 

o.ooo 

0 . ito 

0.114 

0 . lYft 

0.120 

o . aoo 

0 . 144 

o.aas 

0. 191 

o.aao 

0.190 

V . am 

0. 129 

0.900 

0.110 

0.929 

0 oeo 

0.900 

o 002 

0 . 970 

0.070 

O 400 

0 . 007 

0 420 

0 000 

0 410 

0 004 

0 470 

0 002 


0.000 «••- 

e.ooe 


o^/oo 








ORIGINAL PAGE IS 
OK POOR QUALITY 


R/0 
o.ooo 
0 09« 

0 o»o 

0 C7ft 
P 100 

0 . las 
0 lao 
o it% 
o . aop 
0.33B 
O . 2BC 
e . 2Tt 
e SOP 
0 S2I 

0 aco 

o aaa 

0 400 

o aae 

0 «to 
0 a7« 


V»/UP 
O 030 
0 . 03A 
0 039 

o . o«ft 

O OS I 

o . esi 

0.070 
0.07* 
0 . oto 
0 . 0B9 
0 . 003 

o.oss 
0 .oso 
0 o«a 
0 . oso 
0 . OSS 
O . OSS 

o .oas 

O .03S 
e OSS 




'i 

r 


h 







OBIGi^Ai. pfm 
OF POOR QUALITY 


n/o 

0** v«/uo»< 

0 OOO 

0 009730 

0 096 

0 003741 

0 060 

0 002440 

0 0T6 

0 009660 

0 loe 

0 . 004363 

0 136 

0 004667 

0 t6o 

0 013936 

0 176 

0 01361 1 

0 300 

0 013640 

0.33B 

Q 0116*1 

0 360 

0 001706 

0 . 376 

0 007444 

0 . 300 

0 . 006039 

0.336 

0 003066 

0 . 310 

0 003803 

0 376 

0 001410 

0 400 

0 . 001006 

0 436 

0.000633 

0 460 

0 0004 10 

0 476 

0 000266 


0.000 *••>• 
0 OOO 


.010 

u” v*/uo«»a 
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OF POOR QUALITY 


APPENDIX D 

COMPUTER CODE LISTING 
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ORIGiNAL PAGE 
OF POOR quality 


NUMBERED LISTING 


OOOi'^O u 

(ioo?o (,■ **»**»****#*»♦»»♦♦♦♦(■* 4 ******* t*t.*****t*$ 

00000 L » « 

00090 C « THIS F'ftOnRiIM KfmtCFK THE 1 /I IM V-WIRt ItHHiF tiAid IN « 

00100 C » « 

00110 C » NON EUIRI ING H.(IU, JI I Ii I HOI IHF WIRKS « 

00130 C » t 

00130 C » FflLLdU RING'S IrtU. $ 

00140 C * « 

001 SO C « RT « 

00140 C * RRION MCKILLOP * 

00170 C * • * 

00180 C * LAST iteration ON COHIHUR CODE 4-10-03 * 

00190 C « « 

00300 C 

00310 C 
00330 C 

00330 DIMENSION CALU(SO) .CAI I 1 CiO) .CAI I P ( ;.0 > , I nil ( 00 ) r 1 t ( EjO) f ( 3< SO) > 

00340 HADDFC)'?) .SUItF< !iO) • A( I.HE.O' . M I li .1 M'-'IH diO) >ANiil)(00> r 

00350 KFllIf t (PO) rEE'D' i <r;0) - F ,| Ii:'!!' liO » - iMil (S<»> .AM3(S0) »AOl <bO> t 

00340 SAV2<50> .DM ISO • . WVIiO • tIII RHS t '-0 > i"| | fiblHO W HPDHO ( :.() > , VpDHOI 00) » 

00370 SRE YtiTKCvO) rllKPO) rH'HSO) fU 1 :-0 > . l.r.’lillc'f '.O rllUlJOCiO) r VIlUUISO) > 

00380 1! Y C 50 ) . U ( '■,(> ) , ADSls-f; Y f 50 ) 

OO3V0 INTEGFI; DAY. YEAR 

OOjOO RLAL E.DD 

00310 C 

00.330 C ChAPTER 1 

00330 C 5ECM0N 1 - 

00340 C 


00350 C 
00360 C 
00370 C 
00380 
00390 
00400 
00410 C 
00420 C 
00430 C 
00440 C 
004S0 C 
00440 C 
00470 
00400 
<>0490 
00500 
005)0 
00520 
00530 C 
00540 C 
00550 C 
00540 r 
00570 c; 
OOMIO 
00‘n»0 
O0600 
OOf.lO 
00420 
• 00630 
00440 
<)045<.' 
00460 


INFIITING THE NllMIil R tlF PniNTO FRCIM CAL IBRATJ nN.ANOllLATiriN. 
ANU TEST StCTlCiN FXPCWIMFNI 3 .KESF EC T J VELY . 

NMAXolS 

IHAX=3 . 

JHAX=20 


--3FCTIDN 2 


INPiniNG DAI A FROi: CALIDFArUIN EXPERIMENT 

DATA IJ Ai H/2 . 43 . 2 . 76 . 3 . 38 . 3 . 95 . 4 , A 3 . 5 . 1 1 . t,'. 92 . 4 . 72 , 7 , 65 . 
UR .62 .9 , S / , 10 . 7 . 1 1 .8 • . I 3 .05. 1 4 . 4 ' 

READ ( 1 4 . » ) < CALF 1 ' N . . N . 1 . NMAX ) 

RE AD ( 14. 4 > <I'ALE2(N>.N-=1»NMAX) 

URI1E'-6.») <CAI.E1(N> rN-l .NMAX) 

UR I TE f 6 , # > ( CALE 2 < N > . N- 1 f NH A X ) 


SEC II ON 3 

INPIJTII.IH DATA FROM EXPERIMENT IN TF S I '.IFCTiriN 


DATA RDD '<J . 0 . 0 . 025 . 0 • OL, . 0 , 0 ’5 . A 
UO . 25 X5 . O . .3 .<). 325 . 0 . 35 . A . 375 
Rf.AD< 15.4 > (EKJ) » J ) .7MAX) 
KhAD<15r,4) <E2(.|).J 1..IMAX) 

READ ( U. ■ m < ADDE ( J ) , .1.= I , jmAX > 
REAIMIS.*) (fUlDFU)..! I. JMAX) 
REAlH'15-f) <ATTA<.I). I- I. JMAX) 
READt 15.*) < ATI PI J) . I j . JMA.X) 
READ! I5r*) (EHOfi! ( .)' , J 1 .JMAX) 


. I .0. I?5.. 
.0. 4.0. 4 21 


15.0, I 75.0..?.0..-’.:5. 
0.45,0.4/5.0.5/ 
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-y^tfy^rvurjrmtT^r. ■■■ 


S??SSr^qS 


•M'A'^0 UKIth'<^t»> (f|(„().J I.JMi'.X) 

OOAfW WRm<Ap»> (E:M.I).,I |..IMftX> 

00690 URjTr<6tf) <oruir (j> r j-1 • iMi^x ) 

OO/'OO UfiFTr<6.*) (nilRt (J> f.J t» IMftXJ 

OOVIO UKJ 

00730 UWITEFA.*) ! A1 I JK J ) f J-. 1 . JMAX > 

00730 URJ:TE(6.»» <EMVnCi J>rJ=-l»JMAX> 

007-10 r 

0071i0 C KECrrON 1 — 

00760 c 

ooyvo C INriJTINn data FROH ANUUtArinN IXHRIAFNt 

00700 n 

00790 DATA ANOU/l 6 . 1 I ,B7,9 . Ii7t - . , I 1 / 

00800 DATA t IDKE/O . 4U18»0. 47<l-lt 0 . AMT - O . js* J , y . .r.iV// 

00810 DATA tl'DFE70i5B67fO.D461 »0.r,0<l i . 0. ‘Ull’T’ »0 . 

00820 r 

00830 C SECTION fi 

00840 C 

00850 U INPUTINQ HISCEH ANf HUS DAIA 

O0B60 C 

00870 RRAD(16r«) HONTHf DAT r TEAK, XDD -LDD . At i HA .! t. . KAMD , I AMD, EOl ,E02, 

00800 &E01 ,E02 

00*»90 Wf?rTE(6,4) HONTH»riAV.YEAR,XMi,| DD.Al DMA, I T •PAmi, IAMH,E01 ,E02, 

00900 8E01,E02 

00910 C 

00920 C CHAF TF.K I? . 

00930 C 

00940 C *:eciion 1 

00950 C 
00960 C 

0<'970 C THIS STATEMENT TAKE9 CA! H-KATIOH '/fll.TAfiF ONE AND DOES A 

00980 C SErtlNH ORDER LEAST GOUARE!! CURVE FIT Sli|;H THAT! 

00990 t: CALEl*#2=Al+Bl*U**0.5tCl*n 

01000 C 

•MOtO CALL C.‘ I r(NMAy..C7<[-tl,t.UI 1 ,A1 ,D1 ,C1 > 

01070 UK]TE(6. lOOlA! ,B1 ,CI 

01030 JOo FORMAT ( V/70K, 'CALEI SlIHARED ' ,FB . 3, ' + ' ,FH . 3, 2X, 'CALU TO 1/2', 

01040 g +' ,FH.3,2X, 'CALK' ) 

OiOSO C . 

010*0 C SECTION 2 

01070 C 

Oioao r THIS GIATFrtENT TAKES CAI IDTtATrON VDLTAIJt TWO AND DOES A 

OT090 C SECOND ORDEK LEAST SttMAKE!: CURVE F TT Sliai tHAFl 

01 too C PALI :>« • >»A2FDr?*U*»0.ri+C2*lJ 

0 n 1 c 

0117I1 CAIJ. CF 1 MNMAX.CAL ll,cm I 7,A.’,B2,C:2) 

■''1130 _ UIM n '<• . I 10) A2,B7.C:> 

01M<> It',' rtiRMAI I ' ' '30X. 'CAI F7 '3 UHmRE:i ',FH,3.' t ' , FH . 3 ,2X • 'CAUI TO 1/2'* 
01150 S' + ' ,t FI.3,2X» 'CAHJ ' > 

•)I160 C 

01170 t. CHARTER 3 

01180 C 

on°o c no. Lfuir ustno angulation experiment data to compute 

01200 C SENSITIVITY COEFFICIENTS 

01210 r 

0127.0 DO 300 I = 1,IMAX 

.'1,230 ElIi2U< 1 (D1 + '?.0»C1*ANGU( t >*40.5) /< 4.09(A)* ANIHK T ) »Bl»ANnil( I ) »*1 . 5 

OIJ-I.. .’!-tri»'AU(.Ui n*f7)**0.5> 

3125'' t..’li2U( 1 ' tB2F2.0*C2»ANC;iM T)*»0. 5) ' ( 4 . 0* ( A2*AN0Ut 1 >+B2*AMGU( I >»*l .5 

01 260 LFCCfANGIK 1 ) ** 7 >**0 . '.i > 

0JJ7(.‘ URI TF(,.. 1734) FliC'Ui 1 > .F2D.2U( I ) 

0t;;80 1.'34 f MRMATC, /rlOX.F i().6.1(jX.F10.A) 

• 01.290 AMI ( .! > AUniM I ) *E 1 D..>U( I ) 

013OO AM7( 1 ) •O'lU.in |*F21I2U(T) 

0131.. AVI ( J ) I ifiFt < J ) 

OI32(. AV2(1) I J'DFFTT) 


ORlGWAt PA/iS 
OF POOR QUAU«V 


* • 1 \(t 

0! 

Ol.UiO 

.<00 

OI.IAO 

C— 

013 ’0 

c 

01300 

It - 

Oliyo 

c 

OMOO 

l: 

01 A JO 

r. 

01A30 

L 

01 A 30 

C 

01 -MO 
OtAtiO 

11 

01A40 

c 

014 70 
01400 

:?o:? 

01490 

c 

OIDOO 

r 

01510 

c 

01520 

u 

oif.;v.) 

c 

or-. 40 

c 

0151,0 
01 5.', 0 
01570 

:*o.< 

01500 

c 

01590 

c 

01 60A 

<: 

OK., to 

c 

ClA'/O 

r 

01.;, 10 
OlA'IO 
01 olio 
01>.M> 

r' 

oi '0 

»■'••• • 

01 AiiO 

c 

01 A‘/() 

c 

CM VOO 

L 

01 /H) 

i: 

01720 

01 riio 

01 ,’40 

01 ’M) 

Ol/AO 

017,’0 

or/HO 

01 7»0 


OJI'K, 

OlHiO 


Olfi w 
01. '^'40 
OJ USD 
01 H.,ri 
0 1 1 ( '0 
OIHUO 

iO-'O 

■3 1 r ' 0 

Ol .2 
(, t 1 0 

i i :: 

01 y.io 

N f I .< 

0l«4O 

4 

01 /M) 

{• 

01 9 --.0 

f' 

01 ■■■/.. 

C - ’ 

019(10 

f: 


Hh(I ) AMf f 1 )+AM.'< J ) 

noNTiNtir 

CM^fTl R - 

-stcnijN i 


THUi KlflTtMITNI IflKKH Ur- rtNIi mini i\ I INii'iR rUkVI I IT MMI’H I'HAri 
DNeAJ<ia«U 


RtUL C;'F I T ( IHflX , ANBU f DU ■ rt.T • I',< > 

» 3 w 0 . 0 i 32 

URITt(*»202) rt3.B3 

FnRMftT<///20X» 'DM='*F8,3. ' f . Ffl . 1 r 3X r, ' U ' ) 

THIS STDTEMFNI TAKES DU AND Dai: (1 A LINEAR CURVE HT SUCH LHATi 
BtJ=AA + M*U 


stcriON :■ 

CAL U C;’F I T < IMAX r ANOtl . DO . A4 . EH > 

UMTF<e.r303) AAtHA 

FORMA! t /F/ 2 OX 1 'BO=' .t a. ' ♦ ' .F8.3.2Xi 'll' > 

CALI 'll ATING THE Nfl2ZLt EXil VCUICitV 

OOJ = ( ( -H1+(Dl!tT2-A,'.'jiC),'»(Ai -I OILL:') ) T*() . t. 0*0 

UK'’r < , d: H DJ*i2-4.0iC3*( A3't-03*t:') '**0.L)> / <3.0»C2> )**2 

UO=^<LIOL+li02)/2.0 


CHAPTER' A 

DO LOOP TO GENERATE OUTPUT DATA f CAl CUl A T IONS ARE PK'ErUftMEll 
A! LACK R/r. FOSIlinN 
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U'.’f't.KX !P IJ?<..tm'UO 
um‘ 0 (,li= U)IDUO(J) HlSDUOl ,. 1 > )/ 2.0 
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c 
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i: 
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c 

02210 
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r: 
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0.74AO 
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16?0 
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1 

02!, '.-'O 
02!i-'V 
02--. 1-1 
0.24 ' 0 

HO ‘B 
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WKnt( 

FORMArf IHt > 

WRIir snATl HFNin roc tusictl I .-iNI (IDS hafa 

URm:(A.2Sl > MONTH. DAY .yi AR 

t nr<MAT<//,^t,Xf 'ZiAtA tAKI'N' .;!,<• II*. • / ' ,K'. .ir>) 

URnF.(c,,2FJ2> XUD.LOn 

FORMAT! /.;<OY. X/tl«* .F4. 1 .JOX. I /U •' .FA, 1 ) 

alpha. fe 

FORMA I ' / .;i()x. 'ALF'HA .FA. 1 ..‘tiX. FL - ' T. t > 

UR 1 T>: < A • n;<A > famr . t amr . uo . eo i . mfj 

FORMAT !/.2tjX. 'OAhH=' .F4.1 .bX. ' I AMO " .FT., 1 .bX. '110- ' .FS.H. 
85X. *E01«' .F4.2.5X. '£02 = ' ,F4 . 2 ) 


SECTION 

URITf STATEMENT FOR HF:AD1N0H 
URITEiA,2SS> 

FORMA r ( ///.ISX. 'R/D' .7X. 'Ut .'DO' . AX . ' D;>/DO ' . OX . 'U/UO ' .8X . ' UXUO' 
itRX. 'O' UO' f 7X. 'V’/U0'.9X. 'U' W/|l<)**2' ) 

SEC t ION 3 

URin; STATEMENT FOR r.ENEKATED OUTPUT 

WRITE<<..1001 XROD! J) .U1HD0( J) .U2tiU0< J ) .UlUJO( J> .0000(0) . 

XOF 000 ( . ( > . vpimo s ,n . Ri vs r R ( J ) . j= 1 . jmax > 
format '/7.10X.Ftl,3.2X.F0.3.3X.F8.3.;!XrF0.3t4X.F8.3»5X. 

/FD, 3.4*:.r0.3.9’X.Fl().4> 


CHAPTER B 

PLOTTING THE OUTPUT HATA 

URITE<6.e020) 

FORMAT ■! HI) . 

WRITE (A. 8021 ) 

FORMAT (/.30X. 'R/D'.20X. 'U/UO' > 

UR 1 T E ( A . 8022 ) ( RDIi ( ,J ) , UOIIO < J ) . J” 1 . JMAX ) 
F0RHAT(/.27X.F8.3» J5X.F8.3) 

YMA>. -O.S 
XMAX.-3 ,0 

FA( t. F.imPLT(Y.XMAX.YMAX.urillt).Rrin. JMAX) 

w'< n E ' n • ooo<) ) 

I ORM.'l! 'U/UO' ) 

UR VI El. -..0020) 

URITF (4,8023) 

FORMA I ( /.30X, 'R/lf ,20X. 'U/UO' ) 

UE- n E ( /, , 8022 ) ( Y ( J ) , vnUO ( J ) , J= 1 f JMAX ) 
no 1670 J-l.JMAX 
RIlJ'( J)-^ Y( J) 

continue 

CALL POnpL T ( Y . XM AX . YM AX , VIlUO . RIOJ . JMAX ) 
UE rTE«o.800I ) 

> FRMAI ' /'..AX, 'U/UO' > 

UE,nF*-.,.8028) 

EUKHAK ' ,30X. 'R/li ■ .20X. 'U"/UO' ) 

Ul '1 I H ( A • 0022 < I Y ( J ' , OPiniO < J> , J=- 1 . JMAX > 
no i.'. -,! J-l.JMAX 

Rniu. J> Y.J) 

(■oF'ni"Uf 

IV, 1 I, ■ iinfi. T< Y fX.MiTX • VrT,',' .IIFIIIK' ,HIiI> . .IHAX ) 
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ov-'tio Ff<n 

o'2vv(' sufiRnoi iof; criTUxi^Yx^u^FrA^nir.) 

'^*:^ 0 '.o I rm • fi .’inROH mtcs: o srcoNn oRfiCR hirve fit to oiven 
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03700 C 
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03 75 3 r; 
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031100 
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019 to 17 
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(' 1940 
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X4( I ) It'll**;’ 

xm ) 71 1 i*v( I > 

X'.'VI It X’M 1 >*7l I I 
SI I 

sx -sxix< 1 > 

SX2.= 8X2fX;.'( I ) 

SX3 5X33X1( 1 ) 

SX4i.8y4IX4(T) 

BY SYI-Y< I > 

8XY.-SY,y IX Y< I ) 

SX2Y-SX:YfX2l 1 I 1 

URp 1 r,"i5T SfiiiiYPF air.vi ni mf. mnn 

oppi Y PRAMPKs Rim it) <101 Of nil mi'l l t (inori’ifi' t hk 

Aill.ONO c, 

criNriNiiK 

171 --SX;’»gX4-SX.1*SX3 
C2-SX*5X4-SX2*‘5X3 
C3 "SX*i:X 3-SX2*SX? 
oe Y'-s 1 »c I -sx ♦i7:;+sv 2*c3 

B1=EXY»S7 4-S.X;>Y*SX3 
B2-=b'X Y*SX3"SX2Y*SX2 
rtOO:=sY»cmsx*Pi fsx'iyn;’ 

A'=Aflll/IIEr 

D1«SX*HX2Y-5X2*SXY 
AA1 = S 1 »IU -S Y*C2+SX2*ni 
fi" 001/ OKI 

E1«SX*BX2Y-5X2*SXY 
002=SI*(-B2)-SX«El + !iY*r.1 

c-oon/utr 

REriJhN 

end 

SUBRnUTINE C2FIT( rHOXtOMtiU.FONS 1 .iWlD 

rms EUBRoui FNt nors o i.infok miKOf fii hi im ti.-ii.i 
SET nr oni'jU onh cunbi . no f.iiii/.tii'M i > nr m; i .ii". 
UHrPF o ii'i IHE POINT uinrpi <*•.)< y-rNrm rpT > .yni. o !• 
THE SLlIPE 01 THE LINE. 

FHOPTFR I - — 

DIMFNR/riH CIINBT (i.0> lONtilHSO) -X c.iO ) f Y ( '70 ) 

Hxv IS IHE SUM HI X iiMi f; y 

SX rs TMF SUM Of X 

SY IS FHi; SUM OF Y 

S/2 IK THF SUM Ilf X •KIUOKHi 


RXY=0.0 

SX-0.0 

SY”0.0 

SX2»().0 


IIHftPIKR .• 


REPOCMUINO X ONII V 

I'D 100 t l.lM'iX 
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ORlGiNAL f;i5 

OF POOR QUALI1Y 


UNNUMBERED LISTING 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


* * Vr jVi^ iV * A W Vt Vr V!f * * yir yir * yV ^ Vlf Vc yir * jV yif * * *Vlr iV * iHr * * A VnV A ift: * JV jV 3*< * * Vf A * * >V * A 


Ve * 

* THIS PROGRAM REDUCES THE INPUT X-WIRE PROBE DATA IN * 

* , * 

* NON SWIRLING FLOW. IT ISASSUMED THAT THE WIRES * 

ifr 

* FOLLOW KING'S LAW. * 

>V sV 

A by * 

* BRIAN MCKILLOP 

A Vc 

* LAST ITERATION ON COMPUTER CODE 6-10-83 * 

■’t A 
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DIMENSION CALU(SO) ,CALE1 (50) ,CALE2(50) ,RDD(50) ,E1(50) ,E2(50) , 
&ADDE(50) ,SUBE(50) ,ATTA(50) ,ATTB(50) ,EMVDC(50) ,ANGU(50) , 

&E1DFE (50) , E2DFE (SO) , E1D2U (50) , E2D2U (50) , AMl (50) , AM2 (50) , AV 1 (50) , 
4.AV2 (50) , DM (50) . BO (50) , UPRMS (50) , VPRMS (50) , UPDUO (50) , VPDUO (50) , 
&REYSTR(50) ,U1(50) ,U2(50) ,U1DU0(50) ,U2DU0(50) ,UDU0(50) ,VDU0(50) , 
iY(50) ,U(50) ,ABSREY(50) 

INTEGER DAY, YEAR 
REAL LDD 


CHAPTER 1 

SECTION 1 

INPUTING THE NUMBER OF POINTS FROM CALIBRATION, ANGULATION, 
AND TEST SECTION EXPERIMENTS , RESPECT I VELY . 

NMAX=15 

IMAX“5 

JMAX=20 


SECTION 2 


INPUTING DATA FROM CALIBRATION EXPERIMENT 

DATA CALU/2. 43, 2. 76, 3. 38, 3. 95, 4. 63, 5. 11. 5. 92, 6. 72, 7. 65, 

&8. 62, 9. 57, 10. 7, 11. 82, 13. 05, 14. 42, 35*0.0/ 

READ(14,*) (CALEl(N) ,N=1,NMAX) 

READ ( 14,*) (CALE2 (N) , N= 1 , NMaX) 

WRITE (6,*) (CALE1(N),N=>1,NMAX) 

WRITE (6,*) (CALE2(N),N«1,NMAX) 

SECTION 3 

INPUTING DATA FROM EXPERIMENT IN TEST SECTION 

DATA RDD/0. 0,0. 025, 0.05, 0.075, 0.1, 0.125, 0.15, 0.175, 0.2, 0.225, 

&0. 25, 0.275, 0.3, 0.325, 0.35, 0.375, 0.4, 0.425, 0.45, 0.475, 0.5, 29*0.0/ 


112 


nooonoo onnnoor>nor> oonrjo noonr» 


READ05,*) 
READ(15,*) 
READ (15,*) 
READ (15,*) 
READ (15,*) 
READ (15,*) 
READ(15,*) 
WRITE (6,*) 
WRITE (6,*) 
WRITE (6.*) 
WRITE (6.*) 
WRITE (6,*) 
WRITE (6,*) 
WRITE (6,*) 


(E1(J) , J=l, JMAX) 
(E2(J) ,J=1,JMAX) 
(ADDE(J) . J»l, JMAX) 
(SURE (J),J«1, JMAX) 
(ATTA(J) ,J-1,JMAX) 
(ATTB ( J) , J-1 , JMAX) 
(EMVDC(J) , J-1, JMAX) 
(E1(J) , J«l, JMAX) 
(E2(J) , J“l, JMAX) 
(ADDE(J) , J-1, JMAX) 
(SURE ( J) ,J=>1, JMAX) 
(ATTA(J) , J=1,JMAX) 
(ATTB(J) , J"l, JMAX) 
(EMVDC(J) , J=l, JMAX) 
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SECTION A 


INPUTING DATA FROM ANGULATION EXPERIMENT 

DATA ANGU/1A.A2,11,82,9.57,7.65,5.11,A5*0.0/ 

DATA E1DFE/0.61A2, 0.5737,0. 5368, 0.A996,0.A569,A5*0.0/ 
DATA E2DFE/0 . 5592 ,0.5306,0. A938 , 0 . A636 , 0 . AOl 1 , A5*0 . 0/ 

SECTION 5 


INPUTING MISCELLANEOUS DATA 

READ (16,*) MONTH , DAY , TEAR , XDD , LDD , alpha , FE , P AMR , TAMB , EO 1 , E02 , 
&E01,E02 

WR I TE (6 , *) MONTH , DA Y , YEAR , XDD , LDD , ALPHA , FE , PAMB , TAMB , EO 1 , E02 , 
&E01,E02 


CHAPTER 2 
SECTION 1 


THIS STATEMENT TAKES CALIBRATION VOLTAGE ONE AND DOES A 
SECOND ORDER LEAST SQUARES CURVE FIT SUCH THAT: 

CALEl **2-Al+B 1 *U**0 . 5+C 1 *U 

CALL CFI T (NMAX , CALU , CALE 1 , A1 , B 1 , C 1 ) 

WRITE(6,100)A1,B1,C1 

FORMAT (///20X, 'CALEl SQUARED® ’ , F8 . 3 , ' + ' ,F8 .3, 2X, 'CALU TO 1/2', 

&' +' ,F8.3,2X, 'CALU') 

SECTION 2 

THIS STATEMENT TAKES CALIBRATION VOLTAGE TWO AND DOES A 
SECOND ORDER LEAST SQUARES CURVE FIT SUCH THAT: 
CALE2**2»A2+B2*U**0. 5+C2*U 

CALL CFIT(NMAX,CALU,CALE2,A2,B2,C2) 

WRITE (6, 110) A2,B2.C2 
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110 F0RMAT(///20X, 'CALE2 SQUARED^’ .F8.3, ' + ' ,F8. 3, 2X, 'CALU TO 1/2 ' , 

&' +' ,F8.3,2X, 'CALU’) 


C 

C CHAPTER 3 

C 

C DO LOOP USING ANGULATION EXPERIMENT DATA TO COMPUTE 

C SENSITIVITY COEFFICIENTS 

C 


DO 300 I=1,IHAX 

E1D2U(I)'=(B1+2.0*C1*ANGU(I)**0.5)./(4,0*(A1*ANGU(I)+B1*ANGU(I)**1.5 
&+Cl*ANGU (I) *«.’>) **0.5) 

E2D2U (1) <= (B2+2 . 0*C2*ANGU (I) **0 . 5) / 'A . O* (A2*ANGU (I) +B2*ANGU (I) **1 . 5 
&+C2*ANGU ( I) **2) **0.5) 

WRITE(6,1234) E1D2U (I) ,E2D2U(I) 

1234 FORMATC// 10X,F10.6,10X,F10.6) 

AMI (I) “ANGU (I) *E1D2U (I) 

AM2 ( I) =ANGU (I) *E2D2U ( I) 

AV1(I)=E1DFE(I) 

AV2(I)=E2DFE(I) 

DM(I)=AM1(I)+AM2(I) 

BO(l)=AVl(l)+AV2a) 

300 CONTINUE 


C CHAPTER 4 

C 

C SECTION i 

C 

c 

C THIS STATEMENT TAKES DM AND DOES A LINEAR CURVE FIT SUCH THAT; 

C DM-A3+B3*U 

C 

CALL C2FIT(IMAX,ANGU,DM,A3,B3) 

C A3=0.535 

C B3-0.0632 

WRITE (6, 202) A3,B3 

202 F0RMAT(///20X,'DM“',F8.3, ' + ' , F8.3.2X, 'U’) 

C 

C THIS STATEMENT TAKES BO AND DOES A LINEAR CURVE FIT SUCH THAT: 

C B0=A4+B4*U 

C 

C SECTION 2 


C 

CALL C2FIT(1MAX,ANGU,B0,A4,B4) 

WRITE (6, 203) A4.B4 

203 FORMAT (///20X, 'BO=' ,F8. 3, ' +' ,F8.3,2X, 'U') 


C 

C CHAPTER 5 

C 

C CALCULATING THE NOZZLE EXIT VELOCITY 

C 


UOl* ( (-B1+ (B 1 **2-4 . 0*C1 * (Al-E01**2) ) **0 . 5) / (2 . 0*C 1) ) **2 
U02- ( (-B2+ (B2**2-4 . C*C2 * (a2-E02**2) ) **0 . 5) / (2 . 0*C2) ) **2 
U0-(U01+U02)/2.0 
C 

C CHAPTER 6 
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C DO LOOP TO GENERATE OUTPUT DATA .CALCUUTIONS A?.E PREFORMED 

C AT EACH R/D POSITION 
C 

DO 1000 J-1,JHAX 

U1 ( J)“ ( (-B1+ (Bl**2-4. 0*C1* (Al-El ( J) ** 2 ) ) **0. 5) / (2. 0*C1 ) ) **2 

U2 ( J) - ( (~B2+ (B2**2-4 . 0*C2« (A2-E2 ( J) ** 2 ) ) **0 . 5) / (2 . 0*C2 ) ) **2 

U1DU0(J)=U1(J)/U0 

02DU0(J)-U2(J)/U0 

UDUO ( J) * (UlDUO ( J) +U2DU0 ( J) ) /2 . 0 

VDUO ( J) "ABS (UlDUO (J) -U2DU0 ( J) ) /2 . 0 

U(J)»UDU0(J) *U0 

DM(J)-A3+B3*U(J) 

B0(J)»A4+B4*U(J) 

UPRMS ( J) -ADDE ( J) /DM ( J) 

VPRMS (J) -SUBE ( J) /BO O) 

UPDUO (J) -UPRMS (J) *U( J) /UO 
VPDUO ( J) -VPRMS ( J) *U ( J) /UO 

REYSTR ( J) - (EMVDC ( J) *10. 0** ( (ATTA ( J) +ATTB ( J) ) /20. 0) *U ( J) **2/ 

& (U0**2*l 25 . 0*1600 . 0*BO ( J) *DM ( J) ) ) *0 . 00 1 
1000 CONTINUE 

WRITE(6,8113)U0 
8113 FORMAT (/////, 35X.F10. 4) 

C 

C CHAPTER 7 

C 

c sfx:tion 1 

C 

C 

WRITE (6, 250) 

250 FORMAT (IHl) 

C 

C WRITE STATEMENTS FOR MISCELLANEOUS DATA 

C 

WRITE(6,251) month, day, year 

251 FORMAT (//,45X, 'DATA TAKEN* ,2X,I2, V , I2, 7’ .12) 

WRITE (6, 252) XDD.LDD 

252 FORMAT(/,30X, 'X/D-' ,F4,1,30X, *L/D-' ,F4.1) 

WRITE (6, 253) ALPHA. FE 

253 FORMAT (/.30X, 'ALPHA-* ,F4.1,28X, *FE-' ,F4.1) 

WRITE (6, 254) PAMB,TAMB,U0.E01,E02 

254 FORMAT(/,25X, 'PAMB-* ,F4. 1,5X, ’TAMB-’ , F5 . 1 , 5X, 'U0-* ,F8.4, 
C,5X.*E01-' ,F4.2,5X,’E02-’ ,F4.2) 


C 

C SECTION 2 

C 

C WRITE STATEMENT FOR HEADINGS 

C 


WRITE (6, 255) 

255 FORMAT (///,13X, *R/D' ,7X, 'Ul/UO* ,6X, 'U2/U0* ,5X, *U/U0' ,8X, *V/U0 
&9X, 'U"/U0* ,7X, *V*'/U0* ,9X, *U" V"/U0**2’) 

C 

C SECTION 3 

C 
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C WRITE STATEMENT FOR GENERATED OUTPUT 

C 

WRITE(6,1001) (RDD(J) ,U1DU0(J) ,U2DU0U) ,UDU0(J) ,VDU0(J) , 
&UPDU0 ( J) , VPDUO ( J) ,REYSTR ( J) , J-l , JMAX) 

1001 FORMAT(//,10X,F8.3,2X,F8.3,3X,F8.3,2X,F8.3,4X,F8.3,5X, 
&F8.3,4X,F8.3,9X,F10.6) 

C 

Q CHAPTER 8 

C 

C PLOTTING THE OUTPUT DATA 

C 

WRITE (6, 8020) 

8020 FORMAT (IHI) 

WRITE (6, 8021) 

8021 FORMAT (/.30X, 'R/D' ,20X, 'U/UO') 

WRITE (6 , 8022) (RDD ( J) , UDUO ( J) , J=> 1 , JMAX) 

8022 FORMAT(/,27X,F8.3,15X,F8.3) 

YMAX-0.5 

XMAX«1.0 

CALL SUBPLT(Y,XMAX,YMAX, UDUO, RDD, JMAX) 

WRITE (6, 8000) 

8000 FORMAT (/5 AX , ' U/UO ' ) 

WRITE (6. 8020) 

WRITE (6. 8023) 

8023 FORMAT (/,30X, 'R/D' ,20X, 'V/UO') 

WRITE (6, 8022) (Y(J) ,VDU0(J) , J=l, JMAX) 

DO 1670 J-l, JMAX 

RDD(J)-Y(J) 

1670 CONTINUE 

CALL SUBPLT(Y,XMAX,YMAX,VDUO, RDD, JMAX) 

WRITE (6,8001) 

8001 FORMAT (/5AX , ' V/UO ' ) 

WRITE (6, 8020) 

WRITE (6.8028) 

8028 FORMAT (/ , 30X , ' R/D ' , 20X , ’ U"/U0 ' ) 

WRITE (6,8022) (Y(J) ,UPDU0(J) , J-I, JMAX) 

XMAX-0.25 

DO 1673 J-l, JMAX 
RDD(J)-Y(J) 

1673 CONTINUE 

CALL SUBPLT(Y,XMAX,YMAX,UPDUO, RDD, JMAX) 

WRITE(6,800A) 

800A FORMAT (/SAX, 'U"/U0') 

WRITE (6, 8020) 

WRITE (6 ,8029) 

8029 FORMAT (/,30X, 'R/D' ,20X, 'V"/U0') 

WR ITE (6 , 8022) (Y ( J) . VPDUO ( J) , J- 1 , JMAX) 

DO 167A J-l, JMAX 
RDD(J)-Y(J) 

167A CONTINUE 

CALL SUBPLT(Y,XMAX,YMAX, VPDUO, RDD, JMAX) 

WRITE (6, 8005) 

8005 FORMAT (/5AX, 'V"/U0') 

WRITE (6,8020) 
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WRITE (6 8030) 

8030 FORMAT (/, 30X, 'R/D',17X,'U" V'7U0**2') 

WRITE (6, 803 1) (y(J) ,REYSTR(J) , J“l, JMAX) 

8031 FORMAT(/,27X,F8.3,1AX,F10.6) 

DO 1675 J-l,JMAX 
RDD(J)«Y(J) 

1675 CONTINUE 
XMAX-0.02 
DO 1700 J-l.JMAX 
ABSREY ( J) -ABS (REYSTR ( J) ) 

1700 CONTINUE 

CALL SUBPLT(y,XMAX,YMAX, ABSREY, RDD.JMAX) 

WRITE (6, 8006) 

8006 FORMAT (/5 IX, ‘U" V"/U0**2') 

C 

C 

^i((********'ft*ft*******ft*******A***A***/fVf***vr*******yt*ytAff**yf*i(f 
C 

c 

STOP 

END 

SUBROUTINE CFIT(NMAX,U,E,A,B,C) 

C THIS PROGRAM DOES A SECOND ORDER CURVE FIT TO GIVEN 

C HOTWIRE CALIBRATION DATA. 

C THE EQUATION IS OF THE FORM E**2“A+B*SQR (U) +C*U 

C 

C 

DIMENSION E(50) ,U(50) ,Y(50) ,X2(50) ,X3(50) ,XA(50) ,XY(50) , 
*X2Y(50),X(50) 

C 

C CHAPTER 1 

C 

SI=0 

SX“0 

SX2*0 

SX3-0 

SX4-0 

SY-0 

SXY-0 

SX2Y-0 


C 

C CHAPTER 2 

C 

C REPACKAGING X AND Y 

C 


DO 10 I-l.NMAX 

Y(I)-E(I)**2 

X(I)-SQRT(U(D) 

C 

Q CHAPTER 3 

C 

X2(1)-U(I) 

X3(I)-U(I)*X(I) 

XA(I)»U(I)**2 



nnno nnnnooonnnnnno onn 
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XY(I)-X(I)*Y(I) 

X2Y<I)-XY(I)*X(I) 

SI-I 

SX-SX+X(I) 

SX2-SX2+X2(I) 

SX3-SX3+X3(I) 

SX4-SXA+XA(I) 

SY-SY+Y(I) 

SXY-SXY+XY(I) 

SX2Y-SX2Y+X2Y (I) 

USE LEAST SQUARE CURVE FIT METHOD 

APPLY CRAMERS RULE TO SOLVE THE THREE EQUATIONS FOR 


A,B,AND C. 

10 CONTINUE 

C 1 “SX2*SX4-SX3*SX3 
C2“SX*SXA-SX2*SX3 
C3-SX*SX3-SX2*SX2 
DET-SI*C1-SX*C2+SX2*C3 
B1-SXY*SXA-SX2Y*SX3 
B2=SXY*SX3-SX2Y*SX2 
AA0-SY*C1-SX*B1+SX2*B2 
A-AAO/DET 
D1-SX*SX2Y-SX2*SXY 
AA1-SI*B1-SY*C2+SX2*D1 
B-AAl/DET 
E1“SX*SX2Y-SX2*SXY 
AA2-SI* (-B2) -SX*E1+SY*C3 
C-AA2/DET 
RETURN 
END 

SUBROUTINE C2F I T ( IMAX , ANGU , CONST , A, B) 


THIS SUBROUTINE DOES A LINEAR CURVE FIT TO THE DATA 
SET OF ANGU AND CONST. THE EQUATION IS OF THE FORM Y“A+BX 
VIHERE A IS THE POINT WHERE X-O(Y-INTERCEPT) AND B IS 
THE SLOPE OF THE LINE. 


CHAPTER 1 

DIMENSION CONST (50) , ANGU (50) , X (50) , Y (50) 

SXY IS THE SUM OF X TIMES Y 

SX IS THE SUM OF X 

SY IS THE SUM OF Y 

SX2 IS THE SUM OF X SQUARED 

SXY-0.0 

SX-0.0 

SY-0.0 

SX2-0.0 


CHAPTER 2 


REPACKAGING X AND Y 


c 

DO 100 I“1,IMAX 
X(I)=ANGU(I) 

Y(I)=C0NST(I) 

100 CONTINUE 
C 

Q CHAPTER 3 

C 

DO 101 I=1,IMAX 
SXY-SXY+X(I)*Y(I) 

SX“SX+X(I) 

SY=SY+Y(I) 

SX2=SX2+X(I)**2 

101 CONTINUE 
XMEAN=SX/IMAX 
YMEAN“SY/IMAX 
B“1MAX*SXY-SX*SY 
B=B/ (IMAX>'‘SX2-SX**2) 

A“YMEAN-B*XMEAN 

RETURN 

END 

SUBROUTINE SUBPLT ( Y , XMAX , YMAX , XX , Y Y , NLAST) 

C 

C CHAPTER 1 

C 

DIMENSION IM(MOO) ,XX(50) ,YY(50) ,X(50) ,Y(50) 
DO 1200 J=l, NLAST 
X(J) =XX(J) 

Y(J)=YY(J) 

1200 CONTINUE 

WRITE (6, 1500) 

1500 FORMAT (IHI) 

CALL PLOT1(0.6,10,5,20) 

CALL PL0T2(IM,XMAX, 0.0. YMAX, 0.0) 

CALL PL0T3(1H*,X,Y,NLAST,4) 

CALL PLOTAO, 'R/D') 

RETURN 

END 
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SAJffLE ADDRESSED DATA FILES 


U 1 2382A . DECBh . CAL . D 10523 . DATA 


3.487r3.S63»3.631r3*697>3.761r3.817r3.876>3.929»3«988»4.042» 

4*098r4*151f4.210f4*2S9»4*310 

3*502»3*568r3«634»3.700»3.761f3*817r3.873r3*924r3.981r4.032» 

4.088r4*142r4*198f4*247f4.298 


U 1 2382A . DECBM ♦ EXP »D10523.P15.DATA 


4»344r4*344r4»337r4*322r4*318r4«300r4.274»4.244f4.i98r4»125r 

4.061t3*998f3*910»3*807f3*761f3»719r3*639r3.619r3.S95»3«S75 

4»305f4*30Sf4*30S»4.30S>4.296»4.288f4.274r4«242r4.220r4*159» 

4.1i2r4*046r3.973y3«900»3*824r3*768r3.734r3*675f3*631r3*614 

0.0725f0*0742»0.0798»0»0904r0*097sr0«1166f0*139s»0«16i2r 

0 . 1943 f 0 * 2086 r 0 • 2282 1 0 .2288 f 0 . 2384 r 0 . 2431 r 0 . 2396 r 0 . 2447 r 

0.2466*0. 2424r0.2386f 0.2399 

0.0466*0. 0492*0.0510*0. 0589*0. 0685*0.0759*0. 0968*0.1125* 
0.1348*0.1500*0.1629*0.1729*0.1749*0.1725*0.1693*0.1653* 
0.1587*0.1494*0.1386*0.1227 

35.0*35.0*35.0*35.0*36.0*36.0*39.0*39.0*39.0*39.0*39.0*39.0* 

39.0*39.0*39.0*39.0*39.0*39.0*39.0*39.0 

35.0*35.0*35.0*35.0*36.0*36.0*39.0*39.0*39.0*39.0*39.0*39.0* 

39.0*39.0*39.0*39.0*39.0*39.0*39.0*39.0 

254.0*255.0*230.0*277.0*338.0*377.0*520.0*626.0*706.0*734.0* 
715.0*676.0*599.0*505.0*459.0*372.0*324.0*241 .0*186.0*133.0 


U12382A. DECBH. MISC.D10S23.P15. DATA 


5*23*83* 1.5*2.0*90.0*0.0*74.2*27 .0*4.396*4.330*2. 81 *2. 87 


